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The objective of this grant was to obtain accurate potential energy surfaces 
(PES’s) for a number of reactions which are important in the H/N/O combustion 
process. The interest in this at NASA centered around the design of the SCRAM 
jet engine for the National Aerospace Plane (NASP), which was envisioned as an 
air-breathing hydrogen-burning vehicle capable of reaching velocities as large as 
Mach 25. Preli min ary studies indicated that the supersonic flow in the combustor 
region of the scram jet engine required accurate reaction rate data for reactions in 
the H/N/O system, some of which was not readily available from experiment. 

Am ong reactions which are important in the H/N/O system, the following reac- 
tions were selected as being critical and were studied in this grant. (The following 
discussion is organized by class of reaction, while the publications are in chrono- 
logical order.) 

The first class of reactions which were studied are initiation reactions. The 
dominant initiation reaction in H 2 combustion is thought to be the reaction 

J?2 -f- O 2 — * B + HO 2 (1) 

The rate for this reaction has been inferred from the rate of the reverse reaction 

H + HO 2 — ► products (2) 

This analysis is complex because the ractants give rise to both a triplet surface 
which leads to H 2 + O 2 as products and a singlet surface which correlates with 
H 2 O 2 and gives OH as product. The experimental rate is derived from the dif- 
ference in the rate of dissappearance of HO 2 and the rate of appearance of OH. 
There axe considerable uncertainties in the experimental result. Also kinetic mod- 
els used at Lewis and Langley were not in agreement on the rate of this reaction. 
In Ref. 14 the saddle point on the triplet surface for reaction (2) was character- 
ized in sufficient detail to permit computation of the rate constant as a func- 
tion of temperature using transition state theory plus am estimate of tunneling 
through an Eckart barrier. This computed rate is believed to be the most accu- 
rate estimate available for the rate of reaction (2). 

Another initiation reaction which was studied is the reaction 
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OCD) + H 2 -* H 2 0* -» HO + H (3) 

This reaction was of interest because the 1 D excited state of 0 atom reacts with 
H 2 with no barrier, while the 3 P ground state has an « 12.5 kcal/mol barrier. 
Thus, the excited state might be used to initiate combustion in the SCRAM jet. 
This study focused on the long range portion of the potential, in particular on the 
relative barriers to end-on and edge-on insertion of O atom into H 2 . These fea- 
tures of the surface sue relevant to difficulties in reproducing, from calculations, 
the isotopic ratios in the reaction of 0( 1 D) with HD. This work is discussed in 
Ref. 4. 

The most important chain propagation mechsmism in H 2 combustion is the re- 
action 


H + 0 2 ->• H0 2 — *• OH + O (4) 

A detailed study has been made of this reaction. Ref. 3 focused on the minimum 
energy path (MEP) region of this surface. Ref. 8 studied the potential for ex- 
changing the H between the two oxygens via a T-shaped HO 2 saddle point. This 
exchsinge process is important in the reaction of H + O 2 , since the saddle point 
is « 13 kcal/mol below H + O 2 . The study of this system has culminated in the 
last year in the generation of enough computed points to provide a global PES for 
H + O 2 . The H + O 2 and HO 2 portions of the global potential are described in 
Ref. 13, while the OH + 0 portion of the potential is nearly complete and will be 
published later. At the current time the plan is to collaborate with A.F. Wag- 
ner (Argonne National Laboratory) to develop an analytical representation of 
the computed points and carry out dynamics calculations for reaction (4). Fig. 

1 shows a perspective plot of the H + O 2 PES. 

The most important class of combustion reactions from the standpoint of the 
NASP project sire radical recombination reactions, since these reactions result 
in most of the heat release in the combustion process. These processes involve 
recombination of two radicals (A and B) with a third body (C) to remove enough 
energy that the resulting AB fragment remains bound. This process is thought 
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THE H + 0 2 SURFACE 



Fig. 1: Global Potential for H + O 2 . See ref. 13. for a detailed discussion. 


to involve either i) formation of a metastable AB* species which transfers energy 
to C or ii) formation of an intermediate complex AC which subsequently reacts 
with B to yield AB. An example of the latter process, referred to as a chaperone 
mechanism, is illustrated by the series of reactions. 


h + n 2 -> hn 2 

( 5 ) 

H + HN 2 -> H 2 + N 2 

( 6 ) 


The net result of reactions (5) and (6) is recombination of two H atoms with N 2 
as a third body. 
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Calculations described in Refs. 5 and 10 indicate that the HN 2 molecule is 
unstable with respect to H + N 2 by 3-4 kcal/mol but is quasi-bound due to an 
ps 12 kcal/mol barrier to dissociation. The HN 2 well is able to support 6 vibra- 
tional levels and in Ref. 5 the lifetimes for dissociation were estimated using a 
method which utilizes an Eckart barrier to compute one-dimensional tunneling ef- 
fects. In Ref. 10 computed points which define a global potential for reaction (5) 
were reported. The PES for this reaction is shown in Fig. 2. G.S. Schatz (North- 
western) has obtained an analytic representation of the global H + N 2 potential 
and will be carrying out quantum calculations which will determine the rate of 
formation and decay of HN 2 . Calculations for reaction (6) are described in Ref. 

7. The important conclusion for the recombination process is that reaction (6) 
has no barrier and thus should proceed at close to the gas kinetic collision fre- 
quency. Combining this information it should be possible to obtain estimates of 
the rate for the overall recombination process. 

Another species which may be important in H atom recombination is HNO. 

H + NO -► HNO (7) 

H + HNO -+ H 2 + NO ' (8) 

The sequence of reactions (7) and (8) may be important in SCRAM jet simula- 
tion studies due to NO x formation by the Zeldovich mechanism. The lowest 1 A', 
1 A", and 3 A" surfaces for HNO and HON have been characterized in the mini- 
mum energy path regions as discussed in Ref. 9. In a related study, decribed in 
Ref. 11, the rate constant for the reaction, 

NH + O -> N + OH (9) 

which is an important reaction in nitramine combustion, was computed using 
transition state theory with a one-dimensional tunneling correction based on an 
Eckart barrier. While a recommended rate constant expression for reaction (9) 
existed prior to this work, it was only an estimate based on no data and therefore 
the computed rate constant is believed to be much more reliable. Fig. 3 shows 
the computed rate constant for reaction (9). 
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THE H + N 2 SURFACE 



Fig. 2: Global Potential for H + Nj. See ref. 10. for a detailed discussion. 
The last recombination reaction which was considered is 


H + H + H 2 0 - F 2 + H 2 0 (7) 

An analytic function has been developed which consists of the known H 2 and 
H 2 O potentials plus a term which describes the interaction. The interaction term 
is expanded in terms of OH and 00 two body interactions, which were obtained 
from computed H + H 2 O interactions, and a correction term which was fit based 
on H 2 + H 2 O interactions. This potential has an overall root-mean-square error 
of 0.64 mEji. This work is described in Ref. 12. 
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Fig. 3: Computed (transition state theory) rate constant for NH + O — ► N + 

OH. See Ref 11. for a detailed discussion. 
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Appendix 

The appendix contains copies of Ref’s 10-14 which have not yet appeared as jour- 
nal articles. 
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Theoretical Characterization of the Potential Energy Surface 

for 

H + N 2 -► HN 2 

II. Computed Points to Define a Global Potential 

Stephen P. Walch° 

ELORET Institute 
Sunnyvale, Ca. 94087 

Abstract. A previous calculation for H + N 2 [Walch, Duchovic, and Rohlfing, J. 
Chem. Phys, 90, 3230(1989)] focused on the minimum energy path (MEP) region 
of the potential energy surface and on estimates of the lifetime of the HN 2 species. 
In this paper, we report energies computed at geometries selected to permit a global 
representation of the potential energy surface (PES). As in the previous work, the 
calculations were performed using the complete active space self consistent field 
/ externally contracted configuration interaction (CASSCF/CCI) method. The 
surface was characterized using the same basis set as in the previous paper except 
that an improved contraction of the H s basis is used. Calculations with a larger 
basis set were carried out along an approximate MEP obtained with the smaller 
basis set. The new PES exhibits a sharp curvature, which was not present in the 
previous calculations, and has a slightly narrower and smaller barrier to disociation. 
Saddle points for H atom exchange via collinear and T-shaped HN 2 complexes are 
also reported. 

“Mailing Address: NASA Ames Research Center, Moffett Field, CA 94035. 
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I. Introduction 

The HN 2 species has been postulated as an important intermediate in thermal 
De-NO z processes [1-4]. Indirect evidence for the existance of this species has been 
obtained from experimental studies of the reaction: 


NH 2 + NO — > •products (1) 

Studies of reaction (1) have considered three possible product channels [5]: 

N 2 + H 2 O (2) 

N 2 + OH + H (3) 

HN 2 + OH (4) 


While the experimented studies varied widely as to the branching ratio for OH 
production, none of the experiments observed H atoms [5]. This result, along with 
the product channels assumed here, argues for the existence of the HN 2 species. 

Given the possible importance of the HN 2 species in combustion processes, there 
have been several theoretical and experimental studies directed toward estimating 
its lifetime. A recent ab-initio study [6] (hereafter referred to as I) of the minimum 
energy path (MEP) region of the HN 2 surface found the HN 2 species to be unstable 
with respect to H + N 2 by 3.0 keal/mol, but to be quasi-bound due to a 12.2 
keal/mol barrier to dissociation (before zero-point correction). In I the lifetime of 
the HN 2 species was estimated using a method which utilizes an Eckart barrier to 
compute one-dimensional tunneling effects. The lifetime of the lowest vibrational 
level was estimated to be between 8.8 x 10 _n and 5.8 x 10~ 9 sec. This lifetime 
is somewhat longer than the value of 5 x 10~ n obtained theoretically by Curtiss 
et al. [7]. The short lifetime of HN 2 obtained by theory has been supported by 
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the recent experimental work of Selgren et al. [8], who estimate the ground state 
lifetime of HN 2 to be less than 5 X 10” 7 sec. Thus, both theory and experiment 
suggest a short lifetime for HN 2 , which limits its role in combustion processes. This 
is in contrast to combustion models [9] which require lifetimes at least several orders 
of magnitude longer than current estimates. 

Previous theoretical studies of the lifetime of the HN 2 species were based on a 
one- dimensional estimate of tunneling. In order to obtain a more accurate estimate 
of quantum mechanical effects for this reaction a global potential energy surface 
is required. Computations directed toward that goal are reported herein. The 
computational method is discussed in Section II, the results are presented in Section 
III, and the conclusions are given in Section IV. 

II. Computational Details. 

The calculations used the CASSCF/CCI method [10,11] with a selected reference 
list. The details of these calculations are given in I. Two basis sets were used. The 
first basis set was the same as that given in I except for two changes in the H basis 
set. The original H basis used as contraction coefficients the natural orbitals from 
a Cl calculation on H 2 , as described by Almlof and Taylor [12], however using only 
three contracted s type functions leads to an & 1 kcal/mol error in the H atom 
energy, thus we replace the first contracted s type function with the SCF orbital for 
the H atom. In addition, the d function on H was omitted since it was found to have 
only a very small contribution to the energy. The first basis set is [4s3p2dlf/3s*2p] 
and is referred to a s basis set 1, where 3s* indicates the modified contraction of the 
H s functions. The second basis set, which is [5s4p3d2f/4s*3p2d] and is referred to 
as basis set 2, is that given by Almlof and Taylor except for the modified contraction 
of the H s functions. 

The relative positions of the atoms for these calculations are specified in terms of 
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the NN distance (inn), the H to center of mass of N 2 distance ( ih-nn ), and the 
angle (9) between a line connecting H to the center of mass of N 2 and a normal to 
the NN bond at the bond midpoint ( 9 is 0° for T-shaped H-N 2 and 90° for collinear 
H-N 2 ). 

Most of the calculations were carried out in C* symmetry. As discussed below, 
some of the calculations for 0 = 0° and 90° used C 2 „ symmetry. 9 was varied from 
0° through 90° in 10° increments. For each 9 value both ih-nn and inn were 
varied to obtain minimum energy cuts at fixed 9 values. 

The calculations were carried out on the NASA Ames Cray Y-MP/832. These 
calculations used the MOLECULE[13]-SWEDEN[14] system of programs. 

III. Results and Discussion 

The computed energies are given in Table AI of the appendix. In order to aid 
in visualizing the surface, the energy was evaluated along fixed 9 minimum energy 
cuts. For each 6 and ih-nn, t nn was varied and the energy at the minimum 
and the optimal inn are given in Table I. These minimum energy cuts are shown 
graphically in Figs. 1 and 2. 

From Fig. 1 it is evident that for large ih-nn a T-shaped H-N 2 geometry is 
favored. This geometrical arrangement minimizes non-bonded H-N repulsions in 
the long-range repulsive region of the potential. As discussed in ref. 6, at shorter 
ih-nn an NH bond is formed leading to a chemically bonded bent HN 2 species. 
This structure has an L HNN of « 120° (9 « 50°). At ih-nn » 3.5 a 0 , the bent 
HNN structure drops below the T-shaped structure leading to a sharp curvature in 
the reaction path. Formation of HN 2 involves breaking an NN 7r bond as the HN 
bond forms with a resultant barrier. As is evident from Fig. 1, the barrier occurs at 
approximately the same ih—nn as the sharp change in 9. This reflects a crossing 
between the T-shaped structure, which is lowest in the repulsive part of the surface, 
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and the bent HN 2 structure, which is more stable in the chemically bonded region 
of the surface. 

It should be noted here that the PES obtained in I did not show the sharp 
curvature in 9 which is evident in the present PES. This difference appears to be a 
basis set effect. The basis set used in I placed the T-shaped structure too high with 
respect to the bent HN 2 structure and thus the bent structure remains lowest at 
larger th-nn • From Fig. 2 it is evident that the surface is quite flat with respect 
to variations in 9 for r h—NN larger than the ss 4.0 ao. Thus, the sharp curvature 
in 9 observed in the present work may not be especially significant in the dynamics 
of H + N 2 collisions. 

In order to define an approximate MEP, polynomial fits ( six-term quadratic in 
inn and 9) were obtained at each th-nn distance using the three 9 values nearest 
the minimum and three values of inn (9 points). At the HN 2 minimum and en- 
trance channel saddle point ten-term quadratic polynomials in all three coordinates 
were obtained to define the stationary points. It should be noted that this procedure 
defines the stationary points rigorously (within the accuracy of the polynomial fit), 
but the MEP connecting the stationary points is only obtained approximately as 
defined above. The primary problem with the approximate procedure used here to 
define the MEP is that the reaction coordinate is taken as th—nn an d therefore vari- 
ations of the energy along this coordinate are not allowed. In order to indicate that 
this constraint has been imposed we designate the approximate MEP obtained here 
as a constrained energy minimum path (CEM). The CEM is given in Table II. Ad- 
ditional calculations were carried out along the CEM using the [5s4p3d2f/ 4s*3p2d] 
basis set. These results are also given in Table II. In Fig. 3, the energy along the 
CEM with both basis sets is compared with the energy along the CEM of ref. 6 (us- 
ing the basis set of ref. 6). From Fig. 3 it is seen that the present calculations give 
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a narrower barrier and smaller barrier height to HN 2 dissociation. Both of these 
changes result from the improved contraction of the H s basis set which stabilizes 
H + N 2 with respect to HN 2 . 

Table III shows the geometrical parameters of the stationary points on the HN 2 
surface. For the HN 2 minimum and H-N 2 entrance channel saddle point the ge- 
ometries are compared to the results obtained in I. Here it is seen that at the HN 2 
minimum the bond lengths are within 0.01 ao and the bond angle is within 1° of the 
values obtained in I. These differences are probably within the precision of the fits 
used to derive the geometrical parameters. For the H-N 2 minimum, on the other 
hand, the differences are larger. This difference reflects the changes in the large 
th-nn portion of the surface as discussed above. The relative energies are also 
given in Table III. Here it is seen that the best estimate ([5s4p3d2f/4s*3p2d] basis 
set energies at the geometry obtained with the [4s3p2dlf/3s*2p] basis) is that the 
bottom of the HN 2 well is 3.9 kcal/mol above the bottom of the H + N 2 well as 
compared to 3.0 kcal/mol from I. The barrier to dissociation at the same level of 
calculation is 11.3 kcal/mol as compared to 12.2 kcal/mol from I. 

Fig. 4 and Fig. 5 compare the variation of 6 and r nn along the CEM as a 
function of th-nn for the present work and the calculations in I. From Fig. 4 it is 
seen that the $ variation is quite different in the present work, with 6 rising sharply 
from 0° at th-nn about 3.75 a 0 , but for smaller th-nn 0 is similar. On the other 
hand, from Fig. 5 it is seen that the variation in inn is quite similar for the present 
work and ref. I. These results are consistent with the observations made above. 

Figs. 1 and 2 also shows higher energy regions of the surface corresponding to 0 
values near 0° and 90°. Here it is evident that there are saddle points for H atom 
exchange (i.e. motion of an H atom between symmetry equivalent minima on the 
potential energy surface) along cuts with 0 = 0 ° and 90°. For these two choices of 6 
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the actual symmetry is C 2 t> and the barrier (maximum on the surface) prior to the 
saddle point corresponds to a curve crossing ( 2 Aj — * 2 B 2 for 8 = 0° and 2 E + — ► 2 II 
for 8 = 90°). To clarify these regions of the surface, we report energies for both 
electronic states involved in the curve crossing. For 8 = 90 ° the calculations for 
both states were carried out in C 2 i> symmetry leading to an actual curve crossing 
which is evident in Fig. 1 and Table I. For 8 = 90 the 2 B 2 symmetry points were 
computed in C 2 t> symmetry, while the 2 Aj symmetry points were computed in C, 
symmetry. While this region of the surface should also exhibit a curve crossing, this 
region of the surface is plotted in Fig. 1 as if there were an avoided crossing. Since 
this region of the surface is ~ 60 kcal/mol above H + N 2 , this inconsistency should 
not be important. Table III also shows the saddle point geometries and barriers for 
the two exchange saddle points. Both of these saddle points are above the barrier 
height for dissociation of HN 2 to H + N 2 . 

IV. Conclusions. 

Previous studies [6] of the H + N 2 surface have been extended by the use of 
improved basis sets and consideration of larger regions of the surface, with the 
goal of mapping out a globed potential energy surface suitable for use in dynamical 
studies. 

Most of the surface has been characterized with a [4s3p2dlf/3s*2p] basis set which 
differs from the basis set used in ref. 6 only in the contraction of the H s basis. 
This basis provides a more bedanced description of free H atom and bonded H (as 
in HN 2 ). 

In the large r/f-NAT region, the CEM for H atom addition to N 2 obtained with the 
new basis set is different from that obtained in ref. 6. The new CEM shows a sharp 
curvature in the vicinity of the saddle point for H atom addition, with 8 decreasing 
sharply to 0° (T-shaped HN 2 structure) for th-NN greater than 3.5 a<). In the region 
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of the HN 2 minimum the new surface and that of ref. 6 appear to be similar as 
evidenced by essentially identical HN 2 equilibrium geometries. Calculations using a 
[5s4p3d2f/4s*3p2d] basis set along the CEM defined with the [4s3p2dlf/3s*2p] basis 
set place HN 2 3.9 kcal/mol above H + N 2 and predict a barrier to dissociation of 
11.3 kcal/mol (before correction for zero-point energy). These energetics are within 
1 kcal/mol of those predicted in ref. 6. 

The barrier height for dissociation of HN 2 to H + N 2 is slightly smaller and 
the barrier is slightly narrower than in the calculations of ref. 6. Within a one- 
dimensional tunneling model, this result suggests an even shorter lifetime for HN 2 
than had been predicted in ref. 6. 

It should be noted here that one reaction which has not been previously considered 
is: 

H + HN 2 ^H 2 +N 2 ( 5 ) 

Recent studies of the potential energy surface for this reaction [15] indicate no 
barrier for H abstraction. Thus, though the lifetime of HN 2 is short, it is very 
reactive toward H atom, and if formed could effectively scavenge H atoms formed 
in Eq. (3). This could account for the failure to observe H atoms even if both Eq. 
(3) and Eq. (4) are important product channels. 

Saddle points have also been characterized for H atom exchange via T-shaped 
and collinear HN 2 complexes. The barrier heights with the [4s3p2dlf/3s*2p] basis 
set are 46.0 and 29.8 kcal/mol for these two processes, respectively. 
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Table I. H + N 2 energy along fixed 0 cuts®. 


r H-NN 

6 

5.0 

0.0 

4.5 

0.0 

4.0 

0.0 

3.75 

0.0 

3.5 

0.0 

3.25 

0.0 

3.0 

0.0 

2.75 

0.0 

2.25 

0.0 

2.0 

0.0 

1.75 

0.0 

5.0 

10.0 

4.5 

10.0 

4.0 

10.0 

3.75 

10.0 

3.5 

10.0 

3.25 

10.0 

3.0 

10.0 

2.75 

10.0 

2.5 

10.0 

2.25 

10.0 

2.0 

10.0 

1.75 

10.0 

5.0 

20.0 

4.5 

20.0 

4.0 

20.0 


*NN 

Energy 6 

2.1 

- 0.88024 

2.1 

- 0.87748 

2.1 

- 0.87145 

2.1 

- 0.86629 

2.098 

- 0.85890 

2.099 

- 0.84837 

2.102 

- 0.83334 

2.106 

- 0.81194 

2.356 

- 0.79226 

2.369 

- 0.80010 

2.435 

- 0.79515 

2.1 

- 0.88017 

2.1 

- 0.87734 

2.1 

- 0.87127 

2.1 

- 0.86615 

2.099 

- 0.85893 

2.102 

- 0.84885 

2.109 

- 0.83509 

2,126 

- 0.81729 

2.264 

- 0.80302 

2.322 

- 0.81092 

2.351 

- 0.81407 

2.385 

- 0.80219 

2.1 

- 0.87994 

2.1 

- 0.87693 

2.1 

- 0.87073 


6 E(kcal/mol) c 

1.19 

2.92 

6.70 

9.94 

14.58 

21.18 

30.62 

44.04 

56 . 39 * 

51 . 47 * 

54 . 58 * 

1.23 

3.01 

6.81 

10.03 

14.56 

20.88 

29.52 

40.69 

49.64 

44.68 

42.71 

50.16 

1.37 

3.26 

7.15 
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3.75 

20.0 

2.1 

- 0.86573 

10.29 

3.5 

20.0 

2.103 

- 0.85903 

14.50 

3.25 

20.0 

2.111 

- 0.85043 

19.89 

3.0 

20.0 

2.128 

- 0.84046 

26.15 

2.75 

20.0 

2.183 

- 0.83181 

31.58 

2.5 

20.0 

2.260 

- 0.83275 

30.99 

2.25 

20.0 

2.306 

- 0.83593 

28.99 

2.0 

20.0 

2.331 

- 0.82961 

32.96 

1.75 

20.0 

2.374 

- 0.80068 

51.11 

5.0 

30.0 

2.1 

- 0.87950 

1.65 

4.5 

30.0 

2.1 

- 0.87612 

3.77 

4.0 

30.0 

2.1 

- 0.86962 

7.85 

3.5 

30.0 

2.109 

- 0.85898 

14.53 

3.25 

30.0 

2.126 

- 0.85298 

18.29 

3.0 

30.0 

2.162 

- 0.84863 

21.02 

2.75 

30.0 

2.226 

- 0.85129 

19.35 

2.5 

30.0 

2.268 

- 0.85609 

16.34 

2.25 

30.0 

2.297 

- 0.85048 

17.23 

2.0 

30.0 

2.314 

- 0.83173 

31.63 

1.75 

30.0 

2.360 

- 0.77193 

69.15 

5.0 

40.0 

2.1 

- 0.87876 

2.11 

4.5 

40.0 

2.1 

- 0.87474 

4.64 

4.0 

40.0 

2.1 

- 0.86761 

9.11 

3.5 

40.0 

2.119 

- 0.85846 

14.85 

3.25 

40.0 

2.143 

- 0.85564 

16.62 

3.0 

40.0 

2.204 

- 0.85913 

14.43 

2.75 

40.0 

2.224 

- 0.86698 

9.51 

2.5 

40.0 

2.281 

- 0.87015 

7.52 

2.25 

40.0 

2.259 

- 0.85459 

17.28 

2.0 

40.0 

2.263 

- 0.79989 

51.61 
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5.0 

50.0 

2.1 

- 0.87767 

2.80 

4.5 

50.0 

2.1 

- 0.87264 

5.96 

4.0 

50.0 

2.1 

- 0.86425 

11.22 

3.5 

50.0 

2.131 

- 0.85652 

16.07 

3.25 

50.0 

2.177 

- 0.85758 

15.41 

3.0 

50.0 

2.229 

- 0.86608 

10.07 

2.75 

50.0 

2.248 

- 0.87320 

5.60 

2.5 

50.0 

2.237 

- 0.86772 

9.04 

2.25 

50.0 

2.208 

- 0.82346 

36.82 

2.0 

50.0 

2.157 

- 0.71216 

106.66 

5.0 

60.0 

2.1 

- 0.87627 

3.68 

4.5 

60.0 

2.1 

- 0.86980 

7.74 

4.0 

60.0 

2.1 

- 0.85911 

14.45 

3.5 

60.0 

2.141 

- 0.85060 

19.78 

3.25 

60.0 

2.206 

- 0.85444 

17.38 

3.0 

60.0 

2.241 

- 0.86436 

11.15 

2.75 

60.0 

2.242 

- 0.86608 

10.07 

2.5 

60.0 

2.190 

- 0.84011 

26.37 

2.25 

60.0 

2.116 

- 0.75716 

78.42 

5.0 

70.0 

2.1 

- 0.87477 

4.62 

4.5 

70.0 

2.1 

- 0.86655 

9.78 

4.0 

70.0 

2.1 

- 0.85207 

18.86 

3.5 

70.0 

2.147 

- 0.83766 

27.91 

3.25 

70.0 

2.234 

- 0.84331 

24.36 

3.0 

70.0 

2.248 

- 0.85270 

18.47 

2.75 

70.0 

2.209 

- 0.84521 

23.17 

2.5 

70.0 

2.130 

- 0.79685 

53.51 

5.0 

80.0 

2.1 

- 0.87357 

5.37 
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4.5 

80.0 

2.1 

-0.86376 

11.53 

4.0 

80.0 

2.1 

-0.84485 

23.39 

3.5 

80.0 

2.135 

-0.81577 

41.64 

3.25 

80.0 

2.279 

-0.82676 

34.75 

3.0 

80.0 

2.250 

-0.83652 

28.62 

2.75 

80.0 

2.184 

-0.82043 

38.72 

2.5 

80.0 

2.088 

-0.75359 

80.66 

10.0 

90.0 

2.1 

-0.88213 

0.00 

6.05 

90.0 

2.1 

-0.88057 

0.98 

5.05 

90.0 

2.1 

-0.87371 

5.28 

4.55 

90.0 

2.1 

-0.86392 

11.43 

4.05 

90.0 

2.1 

-0.84428 

23.75 

3.75 

90.0 

2.1 

-0.82357 

36.75 

3.5 

90.0 

2.060 

-0.79957 

51. 8 e 

3.5 

90.0 

2.332 

-0.79301 

55.9' 

3.25 

90.0 

2.295 

-0.81616 

41.40 

3.0 

90.0 

2.241 

-0.82552 

35.52 

2.75 

90.0 

2.169 

-0.80590 

47.84 


° Unless otherwise noted energies are in Eh> bond lengths are in ao, and angles are 
in degrees. 

b The energy includes a multireference Davidson’s correction (See Ref. I) and is 
relative to -109.00000 Eh- 
c Energy relative to the H + N 2 asymptote. 

^ These points are for the 2 B 2 state and are computed with C 2 t> symmetry. 
e This point and preceeding points for 6 = 90 0 are for the 2 E + state and are 
computed with C 20 symmetry. 

/ This point and preceeding points for 6 = 90 0 are for the 2 E"*" state and are 
computed with C 2 t> symmetry. 
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Table II. Computed Energies for HN 2 Along the CEM°. 


*H-NN 

ijvn 

0 

Energy (bs l 6 ) 


Energy(bs 2 C ) 


20.5 

2.095 

32.5 

(-.88213) 

0.0 



10.0 

2.1 

0.0 



-109.87596(-. 89360) 

0.0 

5.0 

2.1 

0.0 

(-.88024) 

1.2 

-109.87402(-. 89190) 

1.07 

4.5 

2.1 

0.0 

(-.87748) 

2.9 

-109.87118(-. 88929) 

2.70 

4.0 

2.1 

0.0 

(-.87145) 

6.7 

-109.86501(-.88354) 

6.31 

3.75 

2.1 

0.0 

(-.86629) 

9.9 

-109.85974(-.87860) 

9.41 

3.5 

2.11 

28.4 

(-.85905) 

14.5 

-109.85125(-.87170) 

13.74 

3.35 

2.17 

48.5 

(-.85698) 

15.8 

-109.84416(-. 86936) 

15.21 

3.25 

2.17 

48.9 

(-.85766) 

15.3 

-109.84520(-. 87102) 

14.17 

3.00 

2.23 

52.9 

(-.86636) 

9.9 

-109.85240(-. 88010) 

8.47 

2.75 

2.24 

49.6 

(-.87307) 

5.7 

-109.86002(-.88713) 

4.06 

2.69 

2.25 

47.3 

(-.87321) 

5.6 

-109.86037(-. 88732) 

3.94 

2.5 

2.26 

43.1 

(-.87075) 

7.1 

-109.85655(-. 88656) 

4.42 


0 Unless otherwise noted energies are in Eh, bond lengths are in ao, and singles are 
in degrees. 

b [4s3p2dlf/3s*2p] ANO basis set. Selected reference (C» > 0.05) CASSCF/CCI 
calculations correlating eleven electrons. The values in parenthesis include a mul- 
tireference Davidson’s correction (See Ref. 1) and are relative to -109.00000 Eh- 
The second column of numbers are relative energies in kcal/mol. 
c [5s4p3d2f/4s3p2d] ANO basis set. Selected reference (C* > 0.05) CASSCF/CCI 
calculations correlating eleven electrons. The first column gives the CCI and CCI + 
Davidson’s correction total energies (in Eh), while the second column gives relative 
energies in kcal/mol. 
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Table III. Stationary Points on the H + N 2 Surface 0 . 



H-N 2 

hn 2 

H-N 2 (0=O°) 

H-N 2 (0=9O°) 

rwjv 

2.173 

2.253 

2.390 

2.240 

r H-NN 

3.355 

2.695 

1.955 

3.052 

e 

48.5 

47.3 

0.0 

90.0 

r nh 

2.641 

2.017 



e ,b 

122.7 

115.0 



E(bs 1) 

15.8 

5.6 

51.6 

35.4 

E(bs 2) 

15.2 

3.9 



r i> 

t nn 

2.146 

2.262 



r b 

t nh 

2.753 

2.007 



9 ,b 

118.6 

116.3 



E c 

15.2 

3.0 




° Bond lengths are in ao, angles are in degrees, and energies are in kcal/mol relative 
to the H + N 2 asymptote. 
b L HNN. 

c values from ref. 6. 


Figure Captions. 


Fig. 1. Potential surface for H + N 2 . The figure shows ten fixed 6 minimum energy 
cuts. For each ih-nn, t nn was varied and the minimum energy is shown in the 
figure. 

Fig. 2. Potential surface for H + N 2 . The figure shows the same information as in 
Fig. 1 in the format of a perspective plot. 

Fig. 3. Comparison of energy as a function of tji-nn along the CEM from I and 
the present calculations. For the present calculations results are shown with both 
basis sets. 

Fig. 4. Comparison of 0 as a function of r h-nn for the CEM from I and the present 
work. 

Fig. 5. Comparison of r^vjv as a function of tjj-nn f° r the CEM from I and the 
present work. 
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THE H + N, SURFACE 



/i-B 


NN-Hj 



/6-CL 





Appendix. The appendix contains a table of all the computed CASSCF/CCI ener- 
gies. The energies are in the form CCI(CCI +Q). Note that for the CCI + Q energies 
-109. is not repeated. Thus, for the first point the CCI energy is -109.86500 and 
the CCI+Q energy is -109.88213. Distances are in ao, angles are in degrees, and 
energies are in Eh . 
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Table AI. Computed Energies for HN 2 (C, symmetry) 


r H-NN 

6 

rjviv 

Energy 

20.544 

32.53 

2.095 

-109.86500(-.88213) 

6.587 

37.92 

2.095 

-109.86490(- .88207) 

5.598 

39.33 

2.095 

-109.86378(-.88106) 

5.088 

39.20 

2.096 

-109.86190(- .87934) 

4.608 

40.92 

2.097 

-109.85794(-.87573) 

4.119 

42.24 

2.101 

-109.85037(-.86903) 

3.636 

43.93 

2.118 

-109.83871(-.85970) 

3.446 

44.71 

2.133 

-109.83452(-.85710) 

3.241 

44.37 

2.172 

-109.83273(- .85705) 

3.052 

45.12 

2.205 

-109.83597(- .86164) 

2.872 

46.53 

2.232 

-109.84263(-.86896) 

2.683 

47.38 

2.278 

-109.84633(-.87277) 

2.528 

50.88 

2.255 

-109.83772(-.86807) 

5.00 

0.0 

2.1 

-109.86286(-. 88024) 

4.50 

0.0 

2.1 

-109.85989(-. 87748) 

4.00 

0.0 

2.1 

-109.85347(-. 87145) 

3.75 

0.0 

2.1 

-109.84800(-. 86629) 

3.5 

0.0 

2.0 

-109.83458(-.85234) 

3.5 

0.0 

2.1 

-109.84019(-. 85890) 

3.5 

0.0 

2.2 

-109.83189(-. 85168) 

3.25 

0.0 

2.0 

-109.82335(-.84159) 

3.25 

0.0 

2.1 

-109.82911(-. 84837) 

3.25 

0.0 

2.2 

-109.82098(-. 84139) 
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3.0 

0.0 

2.0 

-109.80740(-.82622) 

3.0 

0.0 

2.1 

-109.81342(-.83334) 

3.0 

0.0 

2.2 

-109.80556(-. 82675) 

2.75 

0.0 

2.0 

-109.78478(-.80428) 

2.75 

0.0 

2.1 

-109.79121(-.81192) 

2.75 

0.0 

2.2 

-109.78378(-.80589) 

2.25 

0.0 

2.2 

-109.76161(-.78499) 6 

2.25 

0.0 

2.4 

-109.76708(-.79170) fc 

2.25 

0.0 

2.6 

-109.74818(-.77464) 6 

2.0 

0.0 

2.2 

-109.76928(-.79222) 6 

2.0 

0.0 

2.4 

-109.77533(-.79983) 6 

2.0 

0.0 

2.6 

-109.75578(-.78536) fc 

1.75 

0.0 

2.2 

-109.75963(-.78234) i> 

1.75 

0.0 

2.4 

-109.76739(-.79486) 6 

1.75 

0.0 

2.6 

-109.74839(-.78885) b 

5.00 

10.0 

2.1 

-109.86278(-.88017) 

4.5 

10.0 

2.1 

-109.85974(-.87734) 

4.00 

10.0 

2.1 

-109.85325(-.87127) 

3.75 

10.0 

2.1 

-109.84779(-.86615) 

3.5 

10.0 

2.0 

-109.83438(-. 85223) 

3.5 

10.0 

2.1 

-109.84011(-.85893) 

3.5 

10.0 

2.2 

-109.83199(-.85191) 

3.5 

10.0 

2.3 

-109.81467(-.83586) 

3.25 

10.0 

2.0 

-109.82340(-. 84176) 

3.25 

10.0 

2.1 

-109.82945(-. 84885) 
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3.25 

10.0 

2.2 

-109.82172(-. 84232) 

3.0 

10.0 

2.0 

-109.80818(-.82717) 

3.0 

10.0 

2.1 

-109.81491(-.83504) 

3.0 

10.0 

2.2 

-109.80809(-. 82951) 

2.75 

10.0 

2.0 

-109.78755(-.80728) 

2.75 

10.0 

2.1 

-109.79593(-.81686) 

2.75 

10.0 

2.2 

-109.79162(-.81383) 

2.5 

10.0 

2.2 

-109.77858(-.80208) 

2.5 

10.0 

2.3 

-109.77631(-.80272) 

2.5 

10.0 

2.4 

-109. 76975(-. 79875) 

2.25 

10.0 

2.2 

-109.77911(-. 80543) 

2.25 

10.0 

2.3 

-109. 78303(-. 81074) 

2.25 

10.0 

2.4 

-109. 77951(-. 80866) 

2.0 

10.0 

2.2 

-109. 77873(-. 80542) 

2.0 

10.0 

2.3 

-109.78458(-.81308) 

2.0 

10.0 

2.4 

-109. 78201(-. 81316) 

1.75 

10.0 

2.2 

-109. 76126(-. 78976) 

1.75 

10.0 

2.3 

-109.76931(-.80008) 

1.75 

10.0 

2.4 

-109. 76984(-. 80213) 

1.75 

10.0 

2.5 

-109. 76474(-. 79837) 

5.00 

20.0 

2.1 

-109.86253(-.87994) 

4.5 

20.0 

2.1 

-109.85926(-.87693) 

4.00 

20.0 

2.1 

-109.85255(-.87073) 

3.75 

20.0 

2.1 

-109.84714(-. 86573) 


20 



3.5 

20.0 

2.0 

-109.83374(-.85187) 

3.5 

20.0 

2.1 

-109.83985(-.85902) 

3.5 

20.0 

2.2 

-109.83224(-.85260) 

3.5 

20.0 

2.3 

-109.81564(-.83737) 

3.25 

20.0 

2.0 

-109.82351(-. 84229) 

3.25 

20.0 

2.1 

-109.83044(-. 85035) 

3.25 

20.0 

2.2 

-109.82393(-.84513) 

3.0 

20.0 

2.0 

-109.81042(-.82999) 

3.0 

20.0 

2.1 

-109.81917(-.83995) 

3.0 

20.0 

2.2 

-109.81514(-.83721) 

2.75 

20.0 

2.1 

-109.80775(- .82938) 

2.75 

20.0 

2.2 

-109.80864(-.83171) 

2.75 

20.0 

2.3 

-109.80209(-.82702) 

2.5 

20.0 

2.2 

-109.8071 0(- .83164) 

2.5 

20.0 

2.3 

-109.80611(-.83227) 

2.5 

20.0 

2.4 

-109.79911(-. 82681) 

2.25 

20.0 

2.2 

-109.80679(-.83203) 

2.25 

20.0 

2.3 

-109.80918(-.83592) 

2.25 

20.0 

2.4 

-109.80413(-.83292) 

2.0 

20.0 

2.2 

-109.79636(-.82300) 

2.0 

20.0 

2.3 

-109.80084(-.82925) 

2.0 

20.0 

2.4 

-109.79820(-.82774) 

1.75 

20.0 

2.2 

-109.76285(-.78997) 

1.75 

20.0 

2.3 

-109.77098(-.79920) 

1.75 

20.0 

2.4 

-109.77114(-.80050) 
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1.75 

20.0 

2.5 

-109.76594(-. 79644) 

5.00 

30.0 

2.1 

-109.86204(-. 87950) 

4.5 

30.0 

2.1 

-109.85834(-. 87612) 

4.00 

30.0 

2.1 

-109.85117(-. 86962) 

3.5 

30.0 

2.0 

-109. 83235(-. 85102) 

3.5 

30.0 

2.1 

-109. 83908(-. 85892) 

3.5 

30.0 

2.2 

-109.83236(-.85353) 

3.5 

30.0 

2.3 

-109.81702(-.83968) 

3.25 

30.0 

2.0 

-109.82331(-.84290) 

3.25 

30.0 

2.1 

-109.83170(-.85256) 

3.25 

30.0 

2.2 

-109.82721(-.84944) 

3.0 

30.0 

2.0 

-109.81354(-.83419) 

3.0 

30.0 

2.1 

-109.82529(-.84720) 

3.0 

30.0 

2.2 

-109.82482(-.84811) 

3.0 

30.0 

2.3 

-109.81668(-.84163) 

2.75 

30.0 

2.1 

-109.82182(-.84478) 

2.75 

30.0 

2.2 

-109.82660(-.85101) 

2.75 

30.0 

2.3 

-109.82314(-.84909) 

2.5 

30.0 

2.1 

-109.82036(-.84411) 

2.5 

30.0 

2.2 

-109.82973(-.85465) 

2.5 

30.0 

2.3 

-109.82964(-.85577) 

2.5 

30.0 

2.4 

-109.82283(-.85065) 

2.25 

30.0 

2.1 

-109.81328(-.83738) 

2.25 

30.0 

2.2 

-109.82531(-.85102) 

2.25 

30.0 

2.3 

-109.82651(-.85468) 
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2.25 

30.0 

2.4 

-109.82141(-.85048) 

2.0 

30.0 

2.2 

-109.80028(-.82697 ) 

2.0 

30.0 

2.3 

-109.80442(-.83166) 

2.0 

30.0 

2.4 

-109.80087(-.82900) 

1.75 

30.0 

2.2 

-109.73752(-.76349) 

1.75 

30.0 

2.3 

-109.74431(-.77111) 

1.75 

30.0 

2.4 

-109.74390(-.77155) 

1.75 

30.0 

2.5 

-109. 73880(-. 76737) 

5.00 

40.0 

2.1 

-109.86124(- .87876) 

4.5 

40.0 

2.1 

-109.85682(-.87474) 

4.00 

40.0 

2.1 

-109.84875(-.86761) 

3.5 

40.0 

2.0 

-109.82968(-.84927) 

3.5 

40.0 

2.1 

-109.83725(-.85822) 

3.5 

40.0 

2.2 

-109.83178(-.85427) 

3.25 

40.0 

2.1 

-109.83258(-.85497) 

3.25 

40.0 

2.2 

-109.83056(-.85441) 

3.25 

40.0 

2.3 

-109.82094(-.84638) 

3.0 

40.0 

2.1 

-109.83138(-. 85494) 

3.0 

40.0 

2.2 

-109.83417(-.85912) 

3.0 

40.0 

2.3 

-109.82908(-.85555) 

2.75 

40.0 

2.1 

-109.83349(-.85775) 

2.75 

40.0 

2.2 

-109.84071(-.86613) 

2.75 

40.0 

2.3 

-109.83901(-.86556) 

2.5 

40.0 

2.1 

-109.83269(-. 85699) 
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2.5 

40.0 

2.2 

-109.84250(-.86782) 

2.5 

40.0 

2.3 

-109.84138(-.87003) 

2.5 

40.0 

2.4 

-109.83351(-.86516) 

2.25 

40.0 

2.2 

-109.82436(-.85343) 

2.25 

40.0 

2.3 

-109.82640(-.85402) 

2.25 

40.0 

2.4 

-109.81962(-.84787) 

2.0 

40.0 

2.2 

-109.77239(-.79853) 

2.0 

40.0 

2.3 

-109. 77264(-. 79941) 

2.0 

40.0 

2.4 

-109.76594(-.79338) 

5.00 

50.0 

2.1 

-109.86009(-.87767) 

4.5 

50.0 

2.1 

-109.85456(-.87264) 

4.00 

50.0 

2.1 

-109.84489(-.86425) 

3.5 

50.0 

2.0 

-109.82496(-.84585) 

3.5 

50.0 

2.1 

-109.83333(-.85592) 

3.5 

50.0 

2.2 

-109.82917(-.85355) 

3.25 

50.0 

2.1 

-109.83116(-.85542) 

3.25 

50.0 

2.2 

-109.83160(-.85738) 

3.25 

50.0 

2.3 

-109.82463(-.85200) 

3.0 

50.0 

2.1 

-109.83408(-.85913) 

3.0 

50.0 

2.2 

-109.83941(-.86573) 

3.0 

50.0 

2.3 

-109.83636(-. 86395) 

2.75 

50.0 

2.1 

-109.83797(-.86292) 

2.75 

50.0 

2.2 

-109.84619(-.87212) 

2.75 

50.0 

2.3 

-109.84466(-.87193) 
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2.5 

50.0 

2.1 

-109.83092(-. 85516) 

2.5 

50.0 

2.2 

-109.83659(-.86682) 

2.5 

50.0 

2.3 

-109.83547(- .86502) 

2.5 

50.0 

2.4 

-109.82699(-. 85629) 

2.25 

50.0 

2.1 

-109.79041(-.81729) 

2.25 

50.0 

2.2 

-109.79623(- .82343) 

2.25 

50.0 

2.3 

-109.79127(- .81891) 

2.0 

50.0 

2.1 

-109.68516(-.71048) 

2.0 

50.0 

2.2 

-109.68518(-. 71118) 

2.0 

50.0 

2.3 

-109.67480(-. 70143) 

5.00 

60.0 

2.1 

-109.85864(-.87627) 

4.5 

60.0 

2.1 

-109.85160(- .86980) 

4.00 

60.0 

2.1 

-109.83931(-.85911) 

3.5 

60.0 

2.0 

-109.81672(-.83878) 

3.5 

60.0 

2.1 

-109.82542(-.84960) 

3.5 

60.0 

2.2 

-109.82221(-.84852) 

3.25 

60.0 

2.1 

-109.82448(-.85027) 

3.25 

60.0 

2.2 

-109.82713(-.85443) 

3.25 

60.0 

2.3 

-109.82238(-.85119) 

3.0 

60.0 

2.1 

-109.82960(- .85539) 

3.0 

60.0 

2.2 

-109.83666(-. 86359) 

3.0 

60.0 

2.3 

-109.83453(-.86281 ) 

2.75 

60.0 

2.1 

-109.83116(-.85612) 

2.75 

60.0 

2.2 

-109.83850(-. 86520) 

2.75 

60.0 

2.3 

-109.83475(-.86445) 
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2.5 

60.0 

2.1 

-109.80891(-. 83542) 

2.5 

60.0 

2.2 

-109.81218(-.84005) 

2.5 

60.0 

2.3 

-109.80424(-.83298) 

2.25 

60.0 

2.0 

-109.72073(-. 74571) 

2.25 

60.0 

2.1 

-109. 73104(-. 75694) 

2.25 

60.0 

2.2 

-109. 72431(-. 75113) 

5.00 

70.0 

2.1 

-109.85713(-.87477) 

4.5 

70.0 

2.1 

-109.84835(-.86655) 

4.00 

70.0 

2.1 

-109.83233(-.85207) 

3.5 

70.0 

2.0 

-109.80363(-.82589) 

3.5 

70.0 

2.1 

-109.81150(-.83647) 

3.5 

70.0 

2.2 

-109.80836(-.83609) 

3.25 

70.0 

2.0 

-109.79430(-.81919) 

3.25 

70.0 

2.1 

-109.80968(-. 83641) 

3.25 

70.0 

2.2 

-109.81468(-.84286) 

3.25 

70.0 

2.3 

-109.81215(-.84167) 

3.0 

70.0 

2.1 

-109.81609(-. 84207) 

3.0 

70.0 

2.2 

-109.82445(-.85158) 

3.0 

70.0 

2.3 

-109.82245(-.85139) 

2.75 

70.0 

2.1 

-109.81248(-.83784) 

2.75 

70.0 

2.2 

-109.81754(-.84516) 

2.75 

70.0 

2.3 

-109.81082(-.84016) 

2.5 

70.0 

2.0 

-109. 75691(-. 78129) 

2.5 

70.0 

2.1 

-109. 76982(-. 79601) 
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2.5 

70.0 

2.2 

-109.76491(-.79239) 

5.00 

80.0 

2.1 

-109.85594(-.87357) 

4.50 

80.0 

2.1 

-109.84566(- .86376) 

4.00 

80.0 

2.1 

-109.82574(-. 84485) 

3.5 

80.0 

2.0 

-109.78759(-.80793) 

3.5 

80.0 

2.1 

-109.79197(-.81525) 

3.5 

80.0 

2.2 

-109.78617(-.81393) 

3.25 

80.0 

2.1 

-109.78609(- .81359) 

3.25 

80.0 

2.2 

-109.79567(-.82460) 

3.25 

80.0 

2.3 

-109.79648(-.82661) 

3.25 

80.0 

2.4 

-109.79009(-.82167) 

3.0 

80.0 

2.1 

-109.79720(-.82374) 

3.0 

80.0 

2.2 

-109.80728(-.83511) 

3.0 

80.0 

2.3 

-109.80537(-.83508) 

2.75 

80.0 

2.1 

-109.78944(-.81562) 

2.75 

80.0 

2.2 

-109.79227(-.82025) 

2.75 

80.0 

2.3 

-109.78164(-.81123) 

2.5 

80.0 

2.0 

-109.72087(-.74564) 

2.5 

80.0 

2.1 

-109.72725(-.75344) 

2.5 

80.0 

2.2 

-109.71316(-.74068) 

21.05 

90.0 

2.1 

- i 09.86483(-.88 i 99 y 

6.05 

90.0 

2.1 

- msessot -. ssos ?) 1 

5.05 

90.0 

2.1 

-109.85613(-.87371) 
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4.55 

90.0 

2.1 

-109.84596(-.86392) b 

4.05 

90.0 

2.1 

-109.82564(-.84428) 6 

3.75 

90.0 

2.1 

-109.80529(-.82357) 1 ’ 

3.5 

90.0 

2.0 

-109.77937(-.79711) 6 

3.5 

90.0 

2.1 

-109.77977(-.79847) 6 

3.5 

90.0 

2.2 

-109.76629(-.78616) 6,c 

3.5 

90.0 

2.2 

-109.75722(-.78543) m 

3.5 

90.0 

2.3 

-109.76266(-.79256) i 

3.5 

90.0 

2.4 

-109.76010(-.79101) 6 

3.25 

90.0 

2.1 

-109.76888(-.79522) fc 

3.25 

90.0 

2.2 

-109.78418(-.81206) 6 

3.25 

90.0 

2.3 

-109.78787(-.81615) 6 

3.25 

90.0 

2.4 

-109.78359(-.81114) 6 

3.0 

90.0 

2.1 

-109.78695(-.81289) b 

3.0 

90.0 

2.2 

-109.79886(-.82445) 6 

3.0 

90.0 

2.3 

-109.79832(-.82330) 6 

2.75 

90.0 

2.1 

-109.77936(-.80255) t ’ 

2.75 

90.0 

2.2 

-109.78223(-.80519) fc 

2.75 

90.0 

2.3 

-109.77041(-.79357) 6 


“ [4s3p2dlf/3s*2p] ANO basis set. Selected reference (Ci > 0.05) CASSCF/CCI 
calculations correlating eleven electrons. 
b Points computed in C 211 symmetry. 

c This point and preceeding points with 6 = 90 0 are 2 E + symmetry. 
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d This point and following points with 0 = 90 ° are 2 II symmetry. 
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Theoretical Characterization of the 
5 H and 3 H Potential Energy Surfaces for 

NH + 0 -» N + OH 


Stephen P. Walch“ 

ELORET Institute 
Sunnyvale, Ca. 94087 

Abstract. The reactant, product, and saddle point regions of the 5 II and 3 II poten- 
tial energy surfaces for the reaction NH + O — ► N + OH have been characterized 
using complete active space self consistent field / externally contracted configu- 
ration interaction (CASSCF/CCI) calculations with large atomic natural orbital 
(ANO) basis sets. The computed barrier heights are 5.6 and 11.7 kcal/mol on the 
5 H and 3 H surfaces, respectively. Transition state theory with an Eckart tunneling 
correction is used to estimate the rate constant on the 5 H surface. 
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I. Introduction 

The lowest three potential energy surfaces (PES’s) of 1 A', 3 A", and 1 A" symmetry 
for H + NO -+ HNO/HON have been discussed elsewhere [1]. In Ref. 1 the barriers 
for H atom addition to NO to give HNO and the regions of the PES’s around the 
HNO /HON eq uili brium geometries were characterized and a schematic of the PES’s 
for HNO /HON was constructed based on ab initio quantum chemical calculations 
[1-4] and qualitative arguments. That schematic of the PES’s for HNO/HON is 
reproduced here as Fig. 1. In this paper, we focus on the PES’s for the reaction of 
NH -I- O. From Fig. 1 it is seen that the reactants NH + O correlate with the 1 A', 
3 A", and 1 A" surfaces of HNO leading to H + NO as one product. 

NH + O — ► HNO -* H + NO ' (1) 

Reaction (1) has been studied by Melius et al. [4] using M0ller-Plesset perturbation 
theory and bond additivity corrections. These authors were mainly interested in 
more complex reactions, which are important in the thermal De-NO* process [5-7], 
but did comment that reaction (1) could proceed with no activation energy on the 
singlet or triplet HNO surfaces. 

From Fig. 1 it is seen that an alternative reaction pathway is the production of 
N + OH. 


NH + O -* N + OH (2) 

Reaction (2) is believed to be important in the combustion of nitramines [8], which 
are of interest as rocket propellants. Benson et al. [9] have recommended the 
expression 
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K = 6.ZE + llT 0S exp(— 8000/ RT)cm 3 mol 1 sec 1 


( 3 ) 


for the rate constant of reaction (2). However, this expression is only an estimate of 
the rate constant and is not based on experimental data. Melius et al. [4] predicted 
a small activation energy for reaction (2) on the quintet surface, where H abstraction 
is the dominant process. In this paper, that abstraction pathway is studied in more 
detail. 

As discussed in section II, reactants and products of reaction (2) may be connected 
by surfaces of 5 H and 3 H symmetry. Of these only the 3 A” component of the 3 II 
surface correlates with the low-lying states of HNO/HON, which have been discussed 
elsewhere [1]. Thus, for bent geometries, the 5 A' and 5 A” components of the 5 13 
surface and the 3 A' component of the 3 H surface do not interact with HNO/HON. 
In addition, for the 3 A" component of 3 H, the collinear geometry corresponds to a 
saddle point. Thus, the 5 H and 3 H surfaces are expected to dominate reaction (2), 
and these surfaces are studied here. 

Qualitative features of the PES’s are discussed in section II, the computational 
method is discussed in Section III, the results axe presented in Section IV, and the 
conclusions are given in Section V. 

II. Qualitative Features. 

The ground state of NH is 3 E~ and the ground state of O is 3 P. From this it 
follows that the ground state reactants give rise to PES’s of 1,3,5 A' and 1,3,5 A" 
symmetry. The ground state of OH is 2 H and the ground state of N is 4 S. Thus, the 
ground state products give rise to PES’s of 3,5 A' and 3,5 A" symmetry. From the 
preceeding discussion, it is seen that PES’s of 3,5 A’ and 3,5 A” symmetry connect 
the reactants and products of reaction (2). For collinear geometries these surfaces 
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correspond to the 3,5 H surfaces considered in this paper. 

The low-lying electronic states of HNO/HON are of 1 A', 1 A", and 3 A" symmetry. 
Thus, only the 3 A" surface of HNO/HON is of the same symmetry as the PES’s for 
reaction (2). As discussed later, the 3 A" component of 3 H is stabilized for highly 
bent geometries (Z NHO less than 135°), but the collinear geometry corresponds to 
a saddle point. From Fig. 1 it is expected that any reactants which “leak” into this 
channel will lead to H + NO, since the barriers to isomerization of HNO to HON 
are larger than the HN bond energy. 

As discussed above, NH + 0 can also lead to surfaces of 1 A' and 1 A" symmetry. 
These surfaces were not studied in this paper, since they can not correlate with the 
products of reaction (2). This channel is expected to lead to H + NO product via 
HNO/HON, as for the 3 A" surface. 

III. Computational Details. 

The calculations are complete active space self consistent field /externally con- 
tracted configuration interaction (CASSCF/CCI) and are basically the same as in 
ref. 1. Two basis sets were used. The smaller basis set is a [4s3p2dlf/3s2pld] atomic 
natural orbital (ANO) basis set as developed by Alml5f and Taylor [10]. This basis 
set is denoted as b.s. 1. and is the same basis set as was used in ref. 1. As discussed 
elsewhere [11-12] b.s. 1 leads to an error of w 1 kcal/mol in the free H atom energy 
and an improved contraction of the H s orbitals has been developed [11]. However, 
for reaction (2) this defect in the basis is not expected to lead to significant errors, 
since free H atom is not involved. As a further check on the accuracy of the results, 
additional calculations were carried out with a larger [5s4p3d2f/4s3p2d] ANO basis 
set, denoted as b.s. 2. 

Most of the calculations were carried out in Cj* symmetry. In the CASSCF 
calculation 10 electrons were correlated (all but the Ols, Nls, and 02s electrons, 
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which leads to a 4<r27r active space. In the subsequent CCI calculations 12 electrons 
were correlated (all but the 01s and Nls electrons). The CAS for eight electrons 
in 3 <t 27 t active orbitals was used as a reference space. 

The saddle points on the 5 II and 3 H surfaces were obtained from a fit to a six-term 
quadratic polynomial in toh and tnh- Bending potentials were obtained for both 
surfaces with r oh and tjvh fixed at the saddle point values. The wavefunctions in 
C, symmetry were equivalent to those for the C2* symmetry calculations, except 
in the case of the 3 II surface only ten electrons were correlated (due to memory 
limitations). 

The calculations were carried out on the NASA Ames Cray Y-MP/832. These 
calculations used the MOLECULE[13]-SWEDEN[14] system of programs. 

IV. Results and Discussion. 

The computed energies for the 5 II and 3 II states are given in Tables I and II, 
respectively. Table III shows the stationary point geometries and energies. In the 
case of the two saddle points the geometry was obtained by fitting to a six-term 
polynomial, while for the reactants and products the geometry was obtained from 
a parabolic fit. 

The computed exoergicity for reaction (2) is 22.9 kcal/mol with b.s 1. Using the 
Do value of 3.37 eV for NH recommended by Bauschlicher et al. [15] in conjunction 
with the experimental [16] Do for OH and experimental [16] u> e values for NH and 
OH, leads to 24.2 kcal/mol as the best estimate of this energy difference. The error 
of 1.3 kcal/mole is twice the error estimate in ref. 15 of ± 0.7 kcal/mol for the Do 
of NH. Calculations with b.s. 2 lead to an exoergicity of 23.1 kcal/mol and reduce 
the error to 1.1 kcal/mol. This basis set differs from that used in ref. 15 in that 
it did not have a g function. This basis set difference and contraction error due to 
the contracted Cl must account for the remainder of the difference between the two 
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calculations. The computed NH and OH bond lengths are « 0.01 a 0 longer than 
experiment. 

The computed barrier heights are 5.6 and 11.7 kcal/mol for the 5 H and 3 H sur- 
faces, respectively, with b.s. 1. The barrier height on the 5 H surface is reduced by 
only 0.07 kcal/mol with b.s. 2. 

Table I also gives energies for bent N-H-0 on the 5 A f and 5 A" surfaces. Table 
IV gives energies for bent N-H-0 on the 3 A' and 3 A" surfaces. In both cases only 
L NHO was varied, with r nh and toh fixed at the saddle point values. In the case 
of the 3 A' and 3 A" surfaces, only ten electrons were correlated (due to memory 
restrictions). This approximation is justified based on the HNO calculations in 
ref. 1, where it was found that the results with 10 and 12 correlated electrons are 
very similar. The bending potentials are also shown in Fig. 2. Here it is seen 
that the 5 A' and 5 A" bending potentials are quite similar. (Note that there is a 
residual splitting for the collinear geometry due to symmetry breaking.) However, 
the bending potentials for the 3 A' and 3 A" PES’s (See Fig. 3.) differ considerably 
for highly bent geometries, with the 3 A" curve dropping for L NHO smaller them 
140°. As discussed earlier, this result is expected, since the 3 A" state correlates 
with HNO/HON for bent geometries. The flat 3 A" bending potential suggests that 
this component of the 3 H surface would lead to H + NO product via HNO*/HON*. 

Table III also gives harmonic frequencies obtained for the saddle point on the 
5 A' and 5 A" surface. These frequencies are based on the polynomial fits discussed 
above and were obtained with the program SURVIB [17]. In order to correct for the 
slight symmetry breaking evident in Fig. 2, the energy at the predicted collinear 
saddle point was computed in C 20 symmetry and the bending curves for the two 
components of the 5 H surface were shifted by the difference between the C 2 * and C, 
sy mm etry calculations. In spite of this correction, there is still a slight (< 4 cm 1 ) 
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difference in the frequencies of the in-plane vibrational modes, which arises from 
different quality least squares fits of the 5 A' and 5 A" energies. 

Transition state theory (TST) calculations were carried out using rigid rotor and 
harmonic oscillator partition functions [18] and an estimate of tunneling through 
an Eckart barrier based on the formalism of Miller [19]. The 5.52 kcal/mol barrier 
obtained with b.s. 2 was used without adjustment. In these calculations the re- 
actants NH and 0 had electronic degeneracy of 3 and spatial degeneracy of 1 and 
3, respectively, while the saddle point N-H-0 species had a spin degeneracy of 5 
(quintet surface) and a spatial degeneracy of 1 for each component ( 5 A* and 5 A” ) 
of the 5 H state. The TST theory calculations were carried out separately for each 
component of the 5 II surface and the rate constants for the two components are 
summed to give the overall rate constant for reaction (2) on the 5 II surface. There 
is also a 5 E surface, but this surface correlates with an excited state of OH and does 
not contribute to the rate for reaction (2). The results of the TST calculations are 
given in Table V. The value of /c, a multiplicative constant which gives the effect 
of tunneling, is also given in Table V. As expected, tunneling is only important 
at lower temperatures: the effect is less than a factor of two for T < 500K. It is 
expected that the 5 H surface will dominate the rate for reaction (2) since the barrier 
for the 3 H surface is twice as large. 

The TST rate constants for the 5 H surface and the rate constant recommended 
by Benson [9] (eqn (3) ) are also plotted as a function of temperature in Fig. 3. 
From Fig. 3 it is seen that Eqn. 3 is a factor of ~ 2.0 below the computed rate 
constant in the high temperature region but is more than two orders of magnitude 
below the computed rate constant at low temperatures. Based on previous work 
[20-21] for similar systems with moderate barriers, the TST theory rate constant is 
expected to be within a factor of 2 of experiment. Because TST does not take into 
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account barrier recrossing effects and because variational effects are not included in 
the present calculation, the TST rate constant should be larger than experiment, 
possibly by as much as the factor of 2 difference between the computed rate constant 
and Eqn. 3 seen at high temperature. However, it is clear that the rate constant at 
low temperature is much larger than predicted by Eqn. 3, and it is believed that 
the TST rate constant is much more reliable and should replace Eqn. 3 as the best 
estimate of the rate constant for reaction (2). 

V. Conclusions. 

The computed barrier heights for the 5 H and 3 H surfaces for 

NH + O -*• N + OH 

are 5.6 and 11.7 kcal/mol, respectively. Thus, the 5 H surface is expected to dom- 
inate the kinetics for this reaction. Transition state theory calculations plus an 
estimate of tunneling based on an Eckart barrier give a rate constant on the 6 H 
surface which is a factor of « 2.0 greater at high temperatures and more than two 
orders of magnitude greater at low temperatures than a rate constant expression 
recommended by Benson (not based on experimental data). While transition state 
theory may slightly overestimate rates ( up to a factor of 2) due to neglect of re- 
crossing effects, it is clear that the computed rate constant is far more reliable than 
the previous recommendation, and should replace it. 

The 3 H surface has also been characterized. The 3 A" component of this surface 
is of the same symmetry as the lower lying HNO/HON surface. This leads to a 
complex bending curve for this component of the 3 H surface and the possibility 
that part of the flux may bleed off into the H + NO product channel. 
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Figure Captions. 


Fig. 1. Schematic of the potential surfaces of HNO/HON (See Ref. 1). Two 
product channels are shown for the reaction of NH + O. H + NO products arise 
from the 1 A', 3 A", and 1 A" surfaces of HNO, while OH + N products arise via an 
H abstraction process predominately on a quintet surface (See the text). 

Fig. 2. Bending potentials for the 5 A' and 5 A" states of N-H-O. 

Fig. 3. Bending potentials for the 3 A ; and 3 A ,f states of N-H-O. 

Fig. 4. Computed rate constant for reaction (2) on the 5 H surface. The rate 
constant is computed separately for the 5 A* and 5 A" components of the 5 H surface 
and the total rate constant is the sum of the separate rate constants. The rate 
expression estimated by Benson is also shown for comparison. 
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Table I. Computed Energies for N-H-0 ( 5 II state) °. 


r NH 

rotf 

6 

Energy 

1.9 

20.0 

180. 

-130. 10584(-. 12822) 

1.994 

20.0 

180. 

-130. 10669(-. 12914) 

2.1 

20.0 

180. 

-130.10403(-. 12653) 

1.9 

2.4 

180. 

-130.08127(-. 10763) 

1.9 

2.6 

180. 

-130.08721(-. 11266) 

1.9 

2.8 

180. 

-130.09191(-. 11639) 

1.9 

3.0 

180. 

-130.09566(-. 11935) 

2.0 

2.4 

180. 

-130.08881(-. 11537) 

2.0 

2.6 

180. 

-130. 09219(-. 11795) 

2.0 

2.8 

180. 

-130.09520(-. 11999) 

2.0 

3.0 

180. 

-130.09792(-. 12187) 

2.1 

2.2 

180. 

-130.09172(-. 11853) 

2.1 

2.4 

180. 

-130.09335(-. 12002) 

2.1 

2.6 

180. 

-130.09413(-. 12016) 

2.1 

2.8 

180. 

-130.09532(-. 12042) 

2.1 

3.0 

180. 

-130.09692(-. 12113) 

2.2 

2.0 

180. 

-130.09485(-. 12106) 

2.2 

2.2 

180. 

-130.09745(-. 12411) 

2.2 

2.4 

180. 

-130.09592(-. 12263) 

2.2 

2.6 

180. 

-130.09409(-. 12034) 

2.2 

2.8 

180. 

-130.09338(-. 11876) 

2.3 

2.0 

180. 

-130.10232(-. 12820) 
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2.3 

2.2 

180. 

-130.10174(-. 12818) 

2.3 

2.4 

180. 

-130.09729(-. 12393) 

2.3 

2.6 

180. 

-130.09290(-. 11926) 

2.3 

2.8 

180. 

-130.09018(-. 11580) 

2.5 

2.0 

180. 

-130.11344(-. 13864) 

2.098 

2.732 

180.* 

-130.09475(-. 12097) 

2.098 

2.732 

170. 

-130.09451(-. 12069) 

2.098 

2.732 

160. 

-130.09374(-. 11984) 

2.098 

2.732 

150. 

-130.09243(-. 11846) 

2.098 

2.732 

140. 

-130.09029(-. 11639) 

2.098 

2.732 

180. c 

-130.09470(*. 12082) 

2.098 

2.732 

170. 

-130.09442(-. 12051) 

2.098 

2.732 

160. 

-130.09358(-.11961) 

2.098 

2.732 

150. 

-130.09211(-. 11806) 

2.098 

2.732 

140. 

-130.08977(-. 11570) 

20.0 

1.7 

180. 

-130.13727(-. 15959) 

20.0 

1.75 

180. 

-130.14087(-. 16327) 

20.0 

1.837 

180. 

-130.14317(-.16572) 

20.0 

1.9 

180. 

-130.14238(-.16503) 


° [4s3p2dlf/3s2pld] ANO basis set. Second order CASSCF/CCI calculations cor- 
relating twelve electrons. 
b Calculation in C, symmetry ( 5 A' state). 
c Calculation in C 4 symmetry( 5 A" state). 
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Table II. Computed Energies for N-H-0 ( 3 II state) °. 


INH 

row 

6 

Energy 

1.9 

2.4 

180. 

-130.07188(-. 09930) 

1.9 

2.6 

180. 

-130. 07911(-. 10506) 

1.9 

2.8 

180. 

-130.08576(-. 11034) 

2.0 

2.2 

180. 

-130.07420(-. 10199) 

2.0 

2.4 

180. 

-130.07816(-. 10581) 

2.0 

2.6 

180. 

-130.08280(-. 10923) 

2.0 

2.8 

180. 

-130.08799(-. 11298) 

2.1 

2.2 

180. 

-130.08147(-. 10898) 

2.1 

2.4 

180. 

-130.08178(-. 10948) 

2.1 

2.6 

180. 

-130. 08349(-. 11039) 

2.1 

2.8 

180. 

-130. 08700(-. 11242) 

2.2 

2.0 

180. 

-130. 08757(-. 11438) 

2.2 

2.2 

180. 

-130.08708(-.11431) 

2.2 

2.4 

180. 

-130.08386(-. 11135) 

2.2 

2.6 

180. 

-130.08240(-. 10966) 

2.2 

2.8 

180. 

-130.08389(-. 10979) 

2.3 

1.8 

180. 

-130.08807(-. 11367) 

2.3 

2.0 

180. 

-130.09522(-. 12150) 

2.3 

2.2 

180. 

-130.09167(-. 11856) 

2.3 

2.4 

180. 

-130.08518(-. 11240) 

2.3 

2.6 

180. 

-130.08045(-. 10783) 

2.3 

2.8 

180. 

-130.07953(-. 10591) 


° [4s3p2dlf/3s2pld] ANO basis set. Second order CASSCF/CCI calculations cor- 
relating twelve electrons. 


15 




Table III. Computed Stationary Point Properties for the 5 II and 3 I3 surfaces of NH 

+ 0 . 




rotf 

6 E 


frequencies 




b.s. 1 

b.s. 2 


NH + 0 

1.975 

20.0 

0.0 

0.0 

3110 

N-H-0( 5 H) 

5 A' 

2.098 

2.732 

5.59 

5.52 

855., 600.(2), -1175. 

5 A" 

2.098 

2.732 

5.59 

5.52 

853., 620.(2), -1179. 

N-H-0( 3 H) 

2.152 

2.571 

11.72 



N + OH 

20.0 

1.847 

-22.9 

-23.1 
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Table IV. Computed Energies for N-H-0 ( 3 II state) Bent Geometries °. 


*NH 

roH 

e 

Energy 

2.152 

2.571 

180. 

-129.99839(-.01616) 

2.152 

2.571 

170. 

-129.99807(-.01536) 

2.152 

2.571 

160. 

-129.99772(-. 01495) 

2.152 

2.571 

150. 

-129.99712(-. 01437) 

2.152 

2.571 

140. 

-129.99626(-. 01360) 

2.152 

2.571 

180. 

-129.99835(-. 01590)' 

2.152 

2.571 

170. 

-129.99817(-. 01566) 

2.152 

2.571 

160. 

-129.99780(-. 01528) 

2.152 

2.571 

150. 

-129.99735(-. 01496) 

2.152 

2.571 

140. 

-129.99712(-.01508) 

2.152 

2.571 

130. 

-129.99790(-. 01663) 


° [4s3p2dlf/3s2pld] ANO basis set. Second order CASSCF/CCI calculations cor- 
relating ten electrons. 
b 3 A' state. 
b 3 A" state. 
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Table V. Computed Transition State Theory Rate Constant for Reaction (1) on the 
5 H Surface. 


T 

K 

5 A' 

rate 

300.00 

3.5326 

0.1265E-13 

500.00 

1.5864 

0.9191E-13 

750.00 

1.2364 

0.3675E-12 

1000.00 

1.1313 

0.9007E-12 

1250.00 

1.0848 

0.1728E-11 

1500.00 

1.0596 

0.2862E-11 

2000.00 

1.0344 

0.6050E-11 

2500.00 

1.0242 

0.1042E-10 

3000.00 

1.0153 

0.1581E-10 


k rate total 

5 A" 


3.5778 

0.1165E-13 

0.2430E-13 

1.5927 

0.8566E-13 

0.1776E-12 

1.2385 

0.3440E-12 

0.7115E-12 

1.1324 

0.8442E-12 

0.1745E-11 

1.0855 

0.1620E-11 

0.3348E-11 

1.0600 

0.2685E-11 

0.5547E-11 

1.0347 

0.5676E-11 

0.1173E-10 

1.0245 

0.9774E-11 

0.2019E-10 

1.0153 

0.1483E-10 

0.3064E-10 


18 


Draft date June 26, 1990 

To be submitted to J. Chem. Phys 


A Potential Energy Surface For H 2 -f H 2 0: 
ab initio calculations and analytical representation. 


David W. Schwenke, 

NASA Ames Research Center, MS 230-3, Moffett Field, CA, 94035 

Stephen P. Walch* and Peter R. Taylor* 

ELORET Institute, Palo Alto, CA 94303 


Abstract: 

We have performed extensive ab initio calculations on the ground state po- 
tential energy surface of H 2 -l-H 2 0, using a large contracted Gaussian basis set and 
a high level of correlation treatment. An analytical representation was then ob- 
tained which represents the calculated energies with an overall root-mean-square 
error which is 0.64 mE h . All nine internal degrees of freedom are explicitly included 
in the fit. The analytical representation is also well behaved as the H 2 bond dissoci- 
ates; it thus can be used to study collision-induced dissociation or recombination of 
H 2 . The strategy used to minimize the number of energy calculations is discussed 
as well as other advantages of the present method for determining the analytical 
represent at ion . 


* Mailing address: NASA Ames Research Center, MS 230-3, Moffett Field, CA, 94035 



I. Introduction 


An important series of reactions for r-odeling combustion are the three-body 
recombination processes 

A + B + M — ► AB + M, (1) 

where A and B are either atoms or small molecules and M is the third body. 
For hydrogen combustion, a wide variety of components have been studied 
experimentally, 1 and a very interesting result is that for A=H, B=OH, most nonre- 
active third bodies give similar rates, with the exception of M=H 2 0, which appears 
to cause a much higher rate of recombination. Because of the experimental uncer- 
tainties involved, it is important to verify this observation with theoretical studies, 
and to explain the cause of the differences. A study of this type by necessity re- 
quires comparison to other third bodies and recombining molecules, and the results 
of extensive calculations for A=B=H, M=H and H 2 have recently been reported. 
The next logical system to consider is A=B=H, M=H20, and in this paper we re- 
port the first step in these calculations, namely the determination of an analytical 
representation of the potential energy of the collision partners. 

The potential energy surface (PES) can be obtained from empirical models, 
from fits to various types of experimental information, from ab initio electronic 
structure calculations, or combinations of these methods. However, for processes 
that sample a large amount of the PES the only practical method is to use ab initio 
data. In the dynamics techniques we are planning to use for the H 2 +H 2 O system, 
we will include highly excited H 2 molecules and collision energies corresponding to 
several thousand degrees Kelvin, 2 thus we are primarily interested in relatively high 
energy portions of surface. Points in arbitrary regions of the PES will be required in 
the dynamics calculations, thus it is necessary to construct an interpolating function 
so that the n umb er of ab initio electronic structure calculations can be kept to a 
minimum. The construction of a faithful interpolating function can be an arduous 
task if an inappropriate strategy is used. 

In this paper we report a successful extension to a nonreactive diatom- 
polyatomic system of a method for constructing interpolating functions originally 
applied to a nonreactive atom- diatom system. 3 This method has also been modified 
for use for a diatom-diatom system . 4 The basic strategy used in the present work 
is to introduce a simple trial function which has many of the global properties of 
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the pentatomic system, and then make corrections to this function to make up for 
its local inadequacies. It is expected that this procedure will reduce the complexity 
of the function which is being represented, resulting in a reduction of the number 
of points required and an increase in the accuracy of the result. The analytical 
representation used in the present work has other attractive features: most of the 
numerical parameters are Unear, which makes their optimization straightforward, 
and also the function is designed for efficient computer evaluation. Many of the 
techniques introduced in this paper concerning polyatomics are more general than 
the present apphcation and are thus expected to be valuable in other attempts 
to produce analytical representations of potential energy surfaces involving poly- 
atomics. 

In Sec. II we discuss our electronic structure calculations, then in Sec. Ill we 
give the interpolating function. Our final discussions and conclusions are contained 
in Sec. IV. All quantities are quoted in Hartree atomic units unless otherwise 
noted, thus the unit of energy is Eh = 2625.500 kJ/mol, the unit of length is 
a 0 = 0.5291771 x 10 -10 m, and the unit of mass is m e = 9.109534 x 10 -31 kg. 

II. Electronic Structure Calculations. 

For a one-electron basis we used a (13s8p6d4f/8s6p4d) Gaussian primitive 
basis set, derived from the van Duijneveldt 5 (13s8p/8s) set by adding polariza- 
tion functions as described in Ref. 6. These functions were then contracted 
using the ANO scheme 6 to yield a [4s3p2dlf/3s2pld] basis set. Finally in or- 
der to improve the description of the outermost regions of the charge density, 
the most diffuse s and p primitive functions on O and the most diffuse prim- 
itive s function on H were uncontracted, 7 ' 8 giving rise to our final basis, de- 
noted [4+ls 3+lp 2d lf/3+ls 2p Id], which was used in all calculations of the 
supermolecule. This basis will give more accurate multipole moments than the 
[4s3p2dlf/3s2pld] basis, and thus can be expected to lead to an improved descrip- 
tion of the long-range interactions. 

The energy calculations were designed to give an accurate description of the 
process H + H + H 2 0 — ► H 2 + H 2 0 with the H*0 geometry always near its 
equilibrium configuration. The most important nondynamical electron correlation 
effects were taken into account by means of CASSCF calculations, and the dynami- 
cal correlation effects were estimated using the size-extensive ACPF method. The 
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CASSCF configurations correspond to a single configuration description of H 2 0 and 
a two-electron two active orbital description of H 2 , and the ACPF calculations used 
the same reference space. All calculations were carried out using the MOLECULE- 
SWEDEN 10 program system on either the NAS CRAY Y-MP/832 or the Ames 
ACF CRAY Y-MP/832. 

The computed H 2 0 fragment properties obtained using these methods are 
given in Table I, where they are also compared to experimental measurements 11 _1S . 
Calculations are carried out for four H 2 0 geometries. These correspond to the 
equilibrium (denoted EQ), 19 an approximate turning point for the symmetric stretch 
normal mode (denoted SS), the asymmetric stretch normal mode (denoted AS), and 
the bending normal mode (denoted B). The normal mode analysis was performed 
using the H 2 0 potential of Ref. 17 which is a slight modification of the fit of Ref. 18 
to the ab initio electronic structure calculations of Ref. 19. The modification is to 
adjust the expansion coefficient a 0 02 to reproduce experimental bending vibrational 
energy levels in a variational basis set calculation. This function represents well the 
results from the present calculation — the function given in Ref. 17 predicts energy 
differences between the equilibrium and distorted geometries given in this table to 
be uni formly 0.1 mEh lower than our ab initio calculations. 

Also in Table I, we compare to experimental values for the dipole moment, the 
qu&drupole moment, and the mean polarizability defined as 

a = — + otyy + a zz ^. (2) 

The quadrupole and dipole moments are computed as expectation values, while 
the polarizability is computed as an energy derivative, as in Ref. 8. For the 
quadrupole moment and mean polarizability, the experimental measurements are 
of vibrationally-thermally-averaged results whereas the theoretical results are for a 
specific geometry. Nonetheless, the agreement is quite satisfactory, with differences 
on the order of a few per cent. For the dipole moment, it is possible to estimate 
values for particular geometries from experiment. To do so we follow Ref. 13 and 
write (note that we use different axis labels and directions than from those in Ref. 
13) 

H'=A+nUi+l*l<l2+l*l 2 ql (3) 

and 

M* = At* 93 » (^) 
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where q> is a dimensionless normal mode coordinate. These coordinates were com- 
puted from 

qi= {]f) £ m i /2 ljiAx ’’ (5) 

where X, is an eigenvalue of the mass-weighted Cartesian Hessian matrix, the sum 
over j is over all coordinates and atoms in H 2 0, m, is the mass of the atom associ- 
ated with index j, Axj is the displacement from equilibrium for the atom/ coordinate 
associated with index j, and /,-< is the j th component of the i th normalized eigenvec- 
tor of the mass-weighted Cartesian Hessian matrix. The eigenvectors were phased 
to be consistent with Ref. 13. For the AS geometry, it is necessary to rotate the 
coordinates before applying Eq. (5) in order to properly separate vibration and ro- 
tation. This is done as described in Ref. 20, and the rotation angle is -5.363°. The 
dipole moment components given by Eqs. (3) and (4) are then rotated by 5.363 
to make the comparisons in Table I. We take from Ref. 11, n\ and from 
Ref. 12, and \i\ from Ref. 13. Again the agreement is satisfactory, with differences 
of a few per cent or less. Of particular note is the fact that the accuracy for the 
different geometries is about the same. 

In Table II we summarize the energies of the noninteracting fragments. We 
will need energies for H 2 bond lengths close to, but not exactly equal to points in 
Table II, and these were generated by interpolation. The six points in Table II with 
H 2 0 at the EQ geometry and H 2 bond length (r Hj ) varying from 1.2 to 2.2 were 
fit to a five term expansion in l/r Hj , and energies were determined at the bond 
lengths 1.201, 1.601, 1.801 and 2.001 a 0 - The difference between using the five term 
expansion and a six term expansion is no larger than 0.01 mEh. 

Four sets of calculations for the interacting system were carried out. The first 
consisted of the calculation of H+H 2 0 interactions at two relative orientations with 
the H 2 0 geometry fixed at the equilibrium geometry. This geometry has OH bond 
lengths equal to 1.8111 a 0 , and an HOH angle equal to 104.45°. These results are 
given in Table III. The second set is made up of 9 different relative orientations of 
H 2 and H 2 0, with the H 2 at its equilibrium separation of 1.401a o and the H 2 0 at 
the experimental equilibrium geometry. The results of these calculations are given 
in Table IV. The third set relaxes the restriction that H 2 be at equilibrium, and the 
results for these calculations are given in Table V. The final set relaxes the H 2 0 
equilibrium geometry restrictions, and the results are given in Table VI. Since the 
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interpolating function obtained using the data of Tables I-V represents the data 
Table VI well, we computed no other points. In all, a total of 181 points were 
computed. It should be noted that although we quote only five decimal digits for 
the energy in the tables, when fitting the points, eight digits were used. Since the 
errors in the fit are larger than O.OlmEh, the difference between using eight and five 
digits should not be significant. 

The relative orientations of H 2 and H 2 0 are given by a code of the form a±bc 
where a is z or x, and b and c axe x, y, or z. This code will uniquely specify the 
orientations in the Cartesian coordinate system which has the 0 at the origin and 
the H 2 center of mass on the positive z axis. The first part of the code gives the 
axis the H 2 bond is parallel to, and the last two elements of the four part code give 
the plane the H 2 0 lies in. Furthermore, the bisector of the HOH angle is the b axis, 
and the sign on b indicates the sign of the b axis components of the coordinates of 
the H atoms in H 2 0. For example, z — zx indicates that the H 2 is on the z axis and 
the H 2 0 lies in the zx plane with the H atoms having coordinates (±1.431545, 0, 
-1.109420). 

III. Interpolating function 

The basic idea of Ref. 3 is that although a pairwise-additive potential is not 
quantitatively accurate, it is always physically reasonable for a nonreactive system, 
thus it can be productive to start out with a pairwise additive potential and then 
make corrections to it. In particular, the pairwise additive potential was shifted by 
a small amoun t to make it always positive, then was multiplied by a function which 
never differs greatly from unity. 

In the extension to a diatom-diatom system, 4 the main modification was to ex- 
plicitly consider the long-range forces and to include the complications arising when 
atom exchange can occur. In the present application we will ignore the possibility 
of atom exchange between the H 2 and H 2 0 and write the potential energy as 

V = V’ nt + V Ha + V H *°, (6) 

where V Hj is the asymptotic H 2 fragment potential, V H, ° is the asymptotic H 2 0 
fragment potential, and V int is the interaction potential. This last term is given by 

V int = s [/ sv (V° + «) - * + V 1H ] + (1 - »)V°, (7) 
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where f^ v is the multiplicative correction function for small vibrational displace- 
ments, j is a switching function, V° is a sum of pairwise potentials, e is a constant 
which ensures that f sv will remain close to unity even if V° is small and V is 
large, or vice versa, and V is the long-range part of the potential. We now discuss 
each of the functions that go into V in detail. 

The function V Ha is taken from Ref. 4, and is an accurate representation of 
the best ground state H2 potential curve. It includes radiative, relativistic and 
adiabatic corrections to the very accurate Born-Oppenheimer calculations of Ref. 
21, although of course these corrections are considerably smaller then our overall 
uncertainties. For the function F Hj °, we take the expression from Ref. 17, which 
is a slight modification of the fit of Ref. 18 to the ab initio electronic structure 
calculations of Ref. 19. The choice of the fragment potentials is not important 
for the determination of the interaction potential and other choices may be more 
suitable for a particular application. For example F H * — ► 00 as r Hj — * 0 for the 
function of Ref. 4 and F Ha ° has spurious deep wells for some geometries highly 
distorted from equilibrium for the function of Ref. 17, and these features may cause 
difficulty. Thus we encourage users of this potential to substitute these fragment 
potentials with other choices if required. 

The zeroth-order potential V° is written as the sum of the two H+H 2 0 inter- 
actions, and each “H3O” potential is written as 

yH.0 , + V HO (fi,o), (8) 

where is a non-bonding pair potential between the lone H and atom A, R HA is 
the distance between the lone H and atom A, and the hydrogens in H 2 0 are labeled 
a and b. The pair potentials are represented as 

F HA = 6 HA exp(-c HA R HA ), (9) 

with and constants. The constants in Eqs. (8)-(9) were determined by 
fitting the data in Table III using nonlinear least squares. For these calculations, 
we take the zero of energy to be — 77.34187Eh, which is the value in Table III 
with the H 2 bond length equal to 10 a 0 . In the least squares procedure, all points 
with interaction energies less than 40 mEh were equally weighted, while the other 
points were given a weight of zero. Thus 17 points were used to determine the 
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four parameters. This fit has a root-mean-square (rms) error of 0.40 mE h and the 
parameters are given in Table VII and the fit is compared graphically to the data 

in Fig. 1. 

The points with interaction energies greater than 40 mEh were not included in 
the fitting procedure because a plot of the energies in Table III for the geometry 
z — zx reveals an avoided crossing for an OH distance ~ 2.0 a<). This feature 
corresponds to an interaction with the l 2 Ai state of H 3 O, which may be thought 
of as a 3 s Rydberg state arising from HaO"*". Since this feature occurs for energies 
greater than 40 mEh, it is not expected to be important in the energy range over 
which the present potential is expected to be accurate. We therefore do not attempt 
to reproduce this feature of the potential, and neglect these higher energy points. 

To proceed further, we must define our coordinates in more detail. Two kinds 
of coordinates are used in the analytical representation of the H3+H3O potential. 
The first consists of distances between points, such as atoms or centers of mass. 
The second is more suitable for scattering calculations and consists of angles as well 
as distances. We will assume that the input to the calculation of the analytical 
representation will consist of Cartesian coordinates of the atoms in some arbitrary 
laboratory-fixed coordinate system. From these coordinates, the calculation of dis- 
tances between points is straightforward, while the determination of the second set 
of coordinates is more complicated. 

The second set of coordinates will consist of $, A, T, r, r Hj , 0, r OH<( , r OH(> , 0, 
a, /?, and 7 . The Euler angles $, A, and T specify the relative orientation between 
the laboratory-fixed coordinate system and the body-fixed coordinate system, r is 
the distance between the centers of mass of H 2 and H 2 O, r Hj is the H 2 bond length, 
6 specifies the H 2 orientation in the body-fixed coordinate system, r OH> , r OBt and © 
specify the H 2 0 geometry, and the Euler angles a, (3 and 7 specify the orientation 
of the water molecule in the body-fixed coordinate system. 

We determine these coordinates from the laboratory-fixed Cartesian coordi- 
nates of the atoms via the following procedure. We first shift the origin of the 
laboratory-fixed coordinates to the center of mass of the H 2 0 molecule, and then 
construct the rotation matrix R Bodv which rotates these laboratory-fixed coordi- 
nates to the body-fixed coordinates. That is, we find the matrix which satisfies 

Body A = \ A (10) 

*i / -j 3 

3 
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and 


x LabA _ £ Rf° dV xf° dyA , (11) 

i 

where xf odyA is the I th component of the Cartesian coordinates of point A in the 
body-fixed coordinate system and xf abA is the j th component of the coordinates of 
point A in the laboratory-fixed coordinate system. Now we choose our body-fixed 
z axis to go from the origin to the center of mass of the H 2 , and the body- fixed xz 
plane to contain H 2 . It is then a simple matter to construct orthogonal unit vectors 
along the body-fixed z and x axes with components expressed the laboratory-fixed 
coordinate system. To get the unit vector along the body-fixed y axis, we simply 
take the cross product of these two unit vectors (y = z x x). Then since the x j 
in Eq. (11) are just columns of the unit matrix, we can identify the elements 
of the rotation matrix R Body with the elements of the unit vectors expressed in 
the laboratory-fixed coordinate system. This rotation matrix determines the Euler 
angles <$, A, T, but we never explicitly need to know these angles it is the rotation 
matrix which is of use. Once R Body has been determined, we can rotate the atomic 
coordinates to the body-fixed coordinate system. 

The lengths r and r Hj are easily determined from the body-fixed coordinates 
of the H 2 , and the angle 6 is the angle which rotates the H 2 about its center of mass 
to lie on the z axis. However as in the case with the angles $,A, and T, it is not 
necessary to explicitly determine 0, but rather just the rotation matrix R Hj which 
rotates the H 2 coordinates so that only the z components are nonzero. The nonzero 
entries of this rotation matrix can be found by simple arithmetic operations from 
the Cartesian coordinates and r Hj . 

The internal H 2 0 coordinates are r OHc , the distance from the 0 to hydrogen 
atom a, , the distance from the 0 to hydrogen atom 6, and 0, the HOH 
angl e. These coordinates are easily determined from the Cartesian coordinates of 
the ato ms and are chosen because they are used for the H 2 0 fragment potential 
of Ref. 17. In order to specify the Euler angles a, /?, 7, we define a standard 
H 2 0 orientation. 23 This consists of the 0 atom lying on the positive z axis, the 
center of mass at the origin, the H atoms in the xz plane, and Hf> having x positive. 
Then the Euler angles are those angles which rotate the particular body frame H 2 0 
coordinates to the standard orientation. As before, it is only necessary to know the 
rotation matrix R Hj ° which rotates the particular H 2 0 coordinates to the standard 
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orientation. This can be determined in the same manner as R Bod », described above, 
without the explicit determination of the Euler angles. 

We now turn to the long-range part of the potential. We will include the 
dipole-quadrupole and quadrupole-quadrupole electrostatic interactions, the dipole- 
induced dipole interaction, and the leading dispersion contribution. These consist 
of terms which decay as r -4 through r -6 . To compute these interactions, we took 
the HjO properties from Table I, the H 2 quadrupole moment from Ref. 24 and 
polarizability from Ref. 25, and the dispersion coefficients from Ref. 26. These 
were then used in damped versions of the appropriate formulas from Ref. 27. 

Specifically, for the dipole-quadrupole interaction, we begin with the expres- 
sions 

= (i2) 

ijk 

T 2 ,ijk = —3 [5 rirjrk-r^Tidjk + rjSik + rkSiji^r- 7 , (13) 

where t, j, k refer to x, y or z components of the various tensors, £i,i is an element of 
the dipole moment tensor for H 2 0, & 2 ,jk is an element of the quadrupole moment 
tensor for H 2 , and r; is a component of the vector from the center of mass of 
the H 2 0 to the center of mass of the H 2 . It should be noted that the multipole 
moments and T{ must all have been calculated using parallel coordinate systems for 
these equations to be valid. We will evaluate Eqs. (12) and (13) in the body-fixed 
coordinate system. In this case the only nonzero contribution from the T{ comes 
from i = 3, (*.e. the z component) thus Eq. (13) becomes 

T 2 ,jjfc = — 3[5^i3^j3^fc3 — ^»3 &jk — &j3$ik ~ • (14) 

We now modify Eq. (14) to remove the singularity at r = 0 by changing r -4 to 
( r 4 + d 4 ) -1 , with d a positive damping parameter. Finally, before we can apply 
Eq. (12), it is necessary to obtain the components of the dipole moment and 
quadrupole moment tensors with respect to the body-fixed coordinates. We do this 
by rotating the laboratory-fixed coordinates using the rotation matrices determined 
when calculating the body-fixed coordinates. That is, in Eq. (12), we use the 
moments 

w.i = E4’°w> (15) 

j 
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ana 

(w) 

kl 

where fij is from Table I, and 

0H’ = -8u{- l) tfh8 (l + h*)Q/ 4, ( 17 ) 

with Q the quadrupole moment from Ref. 24. Now because 0?, J is diagonal and 
R h » only rotates about the y axis, several terms in Eq. (12) vanish, and the so the 
dipole-quadrapole interaction simplifies to 

V = -^2/ii,i02,13 - 3/M,3 02,33]/(r 4 +d 4 ). (18) 

For the quadrupole-quadrupole interaction, we proceed in the same manner, 
starting with 

V @& = 1 ( 19 ) 

y ijfci 

=3 35r<r,-r*ri - 5r 2 {rin8j k + rjr x 8 ik 

+ TkTl8ij + T jT k 8 H + TiTkSjl + TiTjSkl) (20) 

+ r 4 (5i/5jfc + 8ij8 ik + ^fcj^ij)] 7 * 


After including damping and taking advantage of the fact that we evaluate this 
expression in the body frame and that certain elements of the H 2 quadrupole are 
always zero, the quadrupole-quadrupole interaction simplifies to 

F® 0 = 0 / r- JS\ f@2,ll(2®l,ll - 501,33) + 02,22(201,22 “ 501,33) 

3(r 5 + a 5 ) l (21) 

+ 1202,33©1,33 - 1602, 1301, is] » 
where 02,ij Is defined by Eq. (16) and 

<=>,.« = £ °- < 22 > 


with 0® jO from Table I. 

For the dipole-induced dipole interaction we begin with the expression 
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(24) 


i^i° — 2,Jn 


where 

T 2 , it = (3nr t - 6 it r 2 )/r-' 

(25) 

and a 2 ,ij is given by 

&ui = £ r ’"' D « a ‘’' 

kl 

(26) 


where af 3 is an element of the H 2 polarizability tensor. In terms of the || and ± 
components of the polarizability, we have oc ^ = a 22 3 = ol j _ 3 and CK 33 3 = o<|| 3 . The 
other components are zero. Simplifying and damping as before, we obtain 


V ind = -^(« 2 ,nAi,i + 62 , 22 ^ 1,2 + 4 «2,33Ai, 3 ~ + d®). (27) 

Finally, for the dispersion interaction we use 

V di *r = -c 6 /{r* +<£*). (28) 

Then the long-range potential is given by 

V LR = + V e& + V ind + V diap . (29) 

To represent the dependence of the multipole moments, polarizabilities, and 
dispersion coefficients on the vibrational coordinates, we proceeded as follows. For 
the H 2 quadrupole moment, we used the expression from Ref. 4 for Q . This function 
was determined by fitting the accurate results of Ref. 24. 

For the H 2 polarizability, we fit the data of Ref. 25 to the form 

etc = (goc + 9icX + 92cX 2 )/{1 + X 2 )> (30) 

where the subscript c refers to || or J_, gic is a parameter, and 

X = exp[h c (r Hj - r c )]. (31) 

The parameters g ic , h c and r c were determined by nonlinear least squares, with the 
exception that g% c was constrained to be nine so that the asymptotic value of a c 
is accurately obtained. The points were unequally weighted in the fit, with most 
points having a weight of one. The five points near the equilibrium bond length 
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had larger weights, and a point neax the maximum and a point past the maximum 
were also more heavily weighted. The weighted rms error for a|| is 5.7 x 10 a 0 
and the weighted rms error for ax is 4.4 x 10 2 aj. A comparison of these fits and 
the data points is given in Fig. 2, and the values of the parameters are given in 
Table VIII. 

For the dependence of the dispersion coefficient on the H 2 bond length, we 
used the procedure of Ref. 4, i.e. we assume 


c« = C 3 '( a « + a 6 r H, ) «p[-(d< 3 r Hj )’]c 6 , (32) 

with the parameters rrig, dg , 6g, and d.Q taken from Ref. 4. The parameters for 
Eq.(32) are given in Table VIII. 

For the multipole moments and mean polarizability of H 2 0, we make the ap- 
proximation 


M(r 


where M is a component of either the dipole moment, the quadrupole moment, or 
the mean polarizability. The parameters in Eq.(33) are determined by fitting the 
data from Table I and axe given in Table IX. 

To determine the dependence of the dispersion coefficient on the H 2 0 geome- 
try, we invoke the Slater-Kirkwood approximation 28 and the geometric mean rule. 
The Slater-Kirkwood approximation gives the H 2 0-H 2 0 dispersion coefficient as a 
constant times a 3 / 2 , while the geometric mean rule gives the H 2 0-H 2 dispersion 
coefficient as the geometric mean between the H 2 0-H 2 0 and H 2 -H 2 coefficients. 
Thus we take the dispersion coefficient to be proportional to a 3 /* and multiply 
the expression of Eq. (32) by the ratio of the 3/4 power of the mean polarizability 
computed from Eq. (33) and the equilibrium value of 9.35 Uq* 

We now turn to the remaining ingredients of • Before we proceed we need 
to specify the switching function s in Eq.(7). This function is designed to correct 
for deficiencies in the potential when the H 2 bond is greatly stretched. Thus we 
write 

3 = 5/(1 + S), (34) 


with 


3 = exp[-a*(r Hj -r*)]. 


(35) 
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Our procedure to determine f sv will be to guess values for the nonlinear parameters 
e, d, o', and r', and then invert Eq. (7) and solve for the multiplicative correction 
function. This function is then fitted, and then the four nonlinear parameters 
defining the numerically determined f sv are adjusted until a satisfactory fit is 
obtained. 

We will expand the multiplicative correction function as follows: 



/ sv 

pqmm 


,(r,r 


Hj > f OH 4 5 r OH b 


, @)i4pg mm ' ($) & T)j 


pqmm* 


(36) 


where ^pqrnm' ^ ^ yet to be determined function and Apqmm 1 is an angular function. 
The prototype angular function is 29 

Apqmm' = Yp- m (0, 0)2^, (a,/?, 7), (37) 

where Yp m is a spherical harmonic, and 'D^ m i is a Wigner rotation matrix element. 3 
These angular functions are orthogonal, but are complex. This implies relations 
between the expansion coefficients, as f sv is real. Since it can be easily shown that 

Apqmm 1 = ( l)”* Apq—m—m'j (38) 

we can restrict the sum in Eq. (36) to m > 0, m' unrestricted, or m = 0, m' > 0, 
and take the real part of the angular functions. Of the real angular functions with 
q < 2, we have chosen to use the nine angular functions given in Table X. It 
should be noted that for convenience, the normalization factors for these functions 
have been omitted. 

An important simplification is that since we are assuming as input Cartesian 
coordinates the angular functions can be evaluated without trigonometric functions. 
This is because these functions are contained in the rotation matrices. In particular, 

cos(3 = R?i°, (39) 

sin 2 /? cos(2 7 ) = (R*i°) 2 - (R%1°)\ (40) 

cos 4 (/?/2) cos 2(a + 7) = (41) 
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and 


(42) 


sin 4 (/3/2) cos 2(a - 7) 



We now turn to the functions We will first fit the data of Table IV, 

i.e. those points with both the H 2 and H 2 0 at their equilibrium geometry. For this 
fit, the are dependent only on the variable r, and we expand them as 

4-' = M* w*. + £ [-« SJ V - r» + • (43) 

*= 1,3 


Following the principles of Hamilton and Light, 31 we set the exponential parameter 
a SR equal to (1 /d SR ) 2 . In the fitting procedure, for a fixed combination of the six 
nonlinear parameters, r SR , d SR , d, e, a*, and r', we determine the linear parameters 
v sv , ■ that minimize the rms error of the multiplicative correction function f sv . 

pqmm t # , 

The points in this fit are weighted by the inverse of the number of points for a given 
relative orientation, e.g. for geometry z — zx the weight was 1/9 and for geometry 
z + zx, the weight was 1 /8. Then the nonlinear parameters are adjusted to minimize 
the rms errors in the fit to the potential subject to the constraint that r SR > 0. 
After several trials it was determined that the damping parameter d would not be 
large enough to provide adequate damping if it was allowed to vary freely, so it was 
constrained to be always greater than four. This became apparent when H 2 was 
allowed to be distorted away from equilibrium. The switching parameters in Eq. 
(35) axe also adjusted, considering data with the H 2 bond stretched. The resulting 
parameters are given in Table XI. The weighted rms error in the fit to the f sv is 
7.2 x 10 -5 , and the unweighted rms error in V int is only 0.13 mE^. A total of 78 
points are included in this fit. The largest relative errors occur in attractive regions 
of the potential, but since the magnitude of the potential is very small there and we 
plan to use this potential for dynamics calculations at high temperatures, we made 
no attempt to improve the fit to these features by weighting the points differently. 

We next consider fitting the data in Table V, which has includes several values 
of the H 2 bond length. Here we will proceed in a multi-step process. First of 
all, we proceed in a similar manner as described above for the data of Table IV, 
and determine functions (r Hj ) and coefficients Vp^ m m , i( r n 2 )• The current 

procedure differs from that used above in that the nonlinear parameters are fixed 
at the values determined above, and the data of Table V is not used directly. Since 
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slightly different values of H 2 bond lengths were used for the different geometries, 
we synthesize a consistent set of data by taking all points in Tables IV and V which 
differ only in the value of the H 2 bond length, and interpolate to a common group 
of bond lengths. These lengths are 1.2, 1.401, and 1.6 ao, and the interpolation was 
carried out by fitting the points to a polynomial in l/r Hj which exactly reproduced 
the input data. This was only carried out for geometries and values of r where there 
were at least three different values of the H 2 bond length. This yields 35 data points 
per r Hj which were used to determine the 27 parameters v^ mli (r Hj ). These fits to 
f sv gave rms errors in V int from 22^Eh to 62^iEh. Then the resulting parameters 
were fit to the form 

~ ■ < 44 ) 

J=0 

This gives a trial set of expansion coefficients, however it does not exactly reproduce 
the fit to the data of Table IV, since only 35 rather than 78 points axe used for the 
1.401 ao data. We correct this by replacing with the of Table XI. 

This is akin to separately fitting the function and its derivatives.®® The parameters 
for Eq. (44) are given in Table XII. The overall rms error for the 156 data points 
in Tables IV and V with H 2 bond length less than 3 oo using these parameters is 
only 0.64 mEh. 

So far we have not discussed the inclusion of the switching function in Eq. 
(7). The need for such a function can be seen from considering orientations with 
the H 2 on the z axis, such as the z + yx orientation shown in Fig. 3. The ab 
initio calculations predict that as r Hj and r increase with the H to 0 distance fixed, 

decre as es approximately linearly. However, this can not continue indefinitely 
because eventually the further H will no longer interact with the fixed H3 O fragment 
and the interaction energy will become constant. This is exactly what the sum of 
pairwise potentials, V°, predicts, while f sv will continue to predict a decrease. 
Thus we switch to the V° potential for r Hj greater than about 2.5 ao. 

We now turn to the data of Table VI, which has the H 2 at its equilibrium bond 
length and the H 2 0 at distorted geometries. In principle, one could continue the 
procedure used to fit the data in Table V and expand the coefficients in 

terms of the H 2 0 internal coordinates, however this would require a vast increase 
in the number of ab initio energies. We thus proceed in the following manner: the 
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potential as it stands has a dependence on the H 2 0 internal coordinates through 7° 
and V LR , and we ask how well this reproduces the ab initio calculations. Thus we 
performed a series of test calculations with displaced H 2 0 geometries and compared 
to the predictions of the fit. The results of the test calculations are given in Table 
VI, and the comparison between the ab initio and predicted values for V ,nt are 
given in Table XIII. Of the twelve entries in this table, eight have differences less 
than 1 mEh (often much less), three have differences between 1 and 2 mEh, and the 
final point, which is high-energy, has a difference of 8 mEh, which is only 6% of the 
interaction energy. Thus although the fit to this data is not as good as the fit to 
the data in Tables IV and V, the degree of agreement is fair, and thus we conclude 
that the analytical representation is satisfactory as it stands for describing small 
H 2 0 distortions with H 2 at equilibrium. 

We also have not yet considered points where both the H 2 and H 2 0 are dis- 
placed from equilibrium, however, since the analytical representation does a good 
job when either molecule is at its equilibrium geometry, and furthermore no param- 
eters are explicitly required to represent the H 2 0 distortions, we will assume that 
the present analytical representation will describe these points accurately as well. 
Therefore we propose that the analytical representation as it stands will provide a 
realistic representation of the H 2 +H 2 0 interaction potential. 


IV. Discussion and Conclusions 

We have produced an analytical representation of the H 2 +H 2 0 potential energy 
surface suitable for dynamics calculations. All nine internal degrees of freedom are 
included, and the overall rms error between the 156 ab initio points in Tables IV and 
V and the function is only 0.64 mE h . This is obtained with a function containing a 
total of 139 parameters, of which 87 are obtained from the 156 points and 52 which 
are obtained from other data. 

We have computed the minimum energy for the H 2 +H 2 0 van der Waals com- 
plex predicted using this potential. The depth of the minimum is quite small, 0.43 
kcal/mol, which is on the same order of the overall rms error in the fit, 0.40 kcal/mol, 
so the reliability of the prediction can be questioned. However, the rms error for 
attractive points is only 0.08 kcal/mol, so it is likely that the prediction of the van 
der Waals complex is reasonable. The geometry of the minimum is asymmetric, 
but very close to a C 2c symmetry structure, with the H 2 pointing toward the H 2 0 
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center of mass. The angle between the H 2 bond and the H 2 0 plane is about 151 , 
with the H ato ms pointing away from each other. The value of r at the minimum is 
5.81 a 0 and the geometries of the H 2 and H 2 0 are close to their equilibrium values. 

The analytical representation of the H 2 +H 2 0 potential energy surface proposed 
here has several desirable features. It represents the ab initio interaction energies 
well, and does so with relatively few parameters. In addition it is well behaved as H 2 
dissociates, so it can be used to study three-body recombination. From a practical 
point of view, it is useful to observe that most of the least squares parameters are 
linear parameters, so it is easy to determine them. Of the remaining parameters, 
the analytical representation is not a sensitive function of their exact values, so it 
is not difficult to obtain reasonable values for them. From a computational point 
of view, the potential should be reasonably efficient to evaluate, for a minimum 
number of special functions are required. In particular, the only trigonometric 
function used is arc cosine, to determine the HOH angle for the H 2 0 fragment 
potential. From a conceptual point of view, a major advantage of the present 
analytical representation is that it was not necessary to explicitly include parameters 
describing small distortions of the H 2 0 molecule. This greatly reduces the number 
of ab initio calculations required to determine the interpolating function. 
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Table I. Properties of HjO fragment. For H 2 O, the 0 is at the origin and the H 
atoms he in the xz plane. All values quoted are in atomic units. The numbers in 
parenthesis are experimental results. 




H 2 0 geometry 



EQ 

B 

SS 

AS 

±XH 

1.431545 

1.561 

1.515 

1.418 

Z H. 

-1.10942 

-1.013 

-1.233 

-0.9053 


-1.10942 

-1.013 

-1.233 

-1.315 

E acpf 

-76.34202 

-76.33813 

-76.33309 

-76.33295 


-0.7504 

-0.6974 

-0.7762 

-0.7502 


(-0.7268°) 

(-0.6719 6 ) 

(-0.7529 6 ) 

(-0.7258 6 ) 


0. 

0. 

0. 

0.0300 


0. 

0. 

0. 

(0.0304 6 ) 

G C 

v zz 

1.874 

2.382 

1.927 

1.827 


(1.96*) 




0 C 

yy 

-1.757 

-1.875 

-1.922 

-1.759 

(-1.86*) 




0 C 

0. 

0. 

0. 

0.3904 

a e 

9.35 

9.79 

10.39 

9.36 


(9.642^) 





° From Ref. 11. 

6 Computed using the data from Refs. 11 - 13. See text for details. 
c Calculated with the origin at the center of mass. 0« is computed by requiring 
that the quadrupole moment tensor be traceless. 

* From Ref. 14. 

* Mean polarizability computed using the (4+ls 3+lp 2+ld/3+ls 2+lp] basis set 
and methods of Ref. 8. 

f From Ref. 15. 
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Table II. Energies for noninteracting H 2 +H 2 O + 77 Eh- The shortest atom-atom 
distance between the fragments is at least 20 ao . 


r Hj 

H20 b 

e acpf 

1.401 

B 

-0.51168 

1.401 

ss 

-0.50670 

1.401 

AS 

-0.50648 

1.2 

EQ 

-0.50556 

1.401 

EQ 

-0.51555 


_a 

r H, 

H 2 0‘ 

jgACPF 

1.6 

EQ 

-0.50982 

1.8 

EQ 

-0.49633 

2.0 

EQ 

-0.47935 

2.2 

EQ 

-0.46128 

10.0 

EQ 

-0.34187 


0 H 2 bond length. 
fc H 2 O geometry. 
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Table III. Energies for H+H+H 2 0 + 77 E h . All energies and distances are in atomic 


uiubd. 

^OH 

gACPF 

R OH 

jjjACPF 

geometry z — zx b 

geometry z + zx 

1.50 

-0.23800 

2.7000 

-0.27793 

1.75 

-0.28452 

3.2000 

-0.31222 

2.00 

-0.29341 

3.7188 

-0.32870 

2.25 

-0.29812 

3.9688 

-0.33305 

2.50 

-0.30628 

4.2188 

-0.33605 

3.00 

-0.32296 

4.4688 

-0.33810 

3.50 

-0.33308 

4.7188 

-0.33950 

4.00 

-0.33816 

5.2188 

-0.34105 

5.00 

- 0.34149 

5.7188 

-0.34168 

10.00 

- 0.34187 

6.2188 

-0.34190 



7.2188 

-0.34196 

° Distance of closest H to 0 in H 2 0. The other H is 20 a 0 from the 0. 


6 See text for geometry code. 
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Table IV. Energies for H 2 +H 2 0 +77 E h with H 2 and H 2 0 at their equilibrium 
geometries. All energies and distances are in atomic units. 


Rox 

fiACPF 

Rox 

gACPF 


geometry z — zx 


geometry z + z 

2.00 

-0.39854 

2.5 

-0.40682 

2.25 

-0.44006 

3.0 

-0.46738 

2.50 

-0.46767 

3.5 

-0.49552 

2.75 

-0.48595 

4.0 

-0.50775 

3.00 

-0.49787 

5.0 

-0.51467 

3.50 

-0.51022 

6.0 

-0.51555 

4.00 

-0.51479 

7.0 

-0.51558 

5.00 

-0.51650 

20.0 

-0.51554 

10.00 

-0.51564 




geometry z + yx 


geometry x — ; 

2.2 

-0.42505 

2.5 

-0.36057 

2.4 

-0.45114 

3.0 

-0.44664 

2.6 

-0.47025 

3.5 

-0.48662 

2.8 

-0.48412 

4.0 

-0.50384 

3.0 

-0.49410 

5.0 

-0.51378 

3.2 

-0.50119 

6.0 

-0.51527 

3.5 

-0.50807 

7.0 

-0.51547 

4.0 

-0.51361 

20.0 

-0.51554 

5.0 

-0.51616 



10.0 

-0.51561 




geometry x — zy 


geometry x + 

2.5 

-0.36970 

2.5 

-0.22643 

3.0 

-0.45060 

3.0 

-0.38570 

3.5 

-0.48848 

3.5 

-0.45984 

4.0 

-0.50480 

4.0 

-0.49283 

5.0 

-0.51408 

5.0 

-0.51267 

6.0 

-0.51539 

5.5 

-0.51486 
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7.0 


-0.51552 

6.0 

-0.51561 

20.0 


-0.51555 

6.5 

-0.51582 




7.0 

-0.51583 




20.0 

-0.51555 


geometry 

x + zy 


geometry x +yx 

2.5 


-0.28418 

2.5 

-0.34296 

3.0 


-0.40844 

3.0 

-0.43577 

3.5 


-0.46952 

3.5 

-0.48047 

4.0 


-0.49730 

4.0 

-0.50060 

5.0 


-0.51388 

5.0 

-0.51307 

5.5 


-0.51556 

6.0 

-0.51521 

6.0 


-0.51604 

7.0 

-0.51552 

7.0 


-0.51601 

20.0 

-0.51555 

20.0 


-0.51555 




geometry 

x + xy 


geometry x + xy 

2.75 


- 0.40140 

4.00 

-0.50137 

3.00 


-0.43867 

4.50 

-0.50980 

3.25 


-0.46443 

5.00 

-0.51334 

3.50 


-0.48190 

10.00 

-0.51555 

a Distance between 0 and nearest H if the H 2 bond points to the 0, otherwise the 


distance between 0 and the H 2 center of mass. 
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Table V. Energies for H 2 +H 2 0 + 77 E h with H 2 0 at its equilibrium geometry and 
H 2 displaced from equilibrium. All energies and distances are in atomic units. 


*ox 


e acpf 

R ox 


gACPF 

geometry 

z — zx 

geometry 

x + xy 

2.00 

1.6 

-0.41306 

2.75 

1.601 

-0.39807 

2.00 

1.8 

-0.41768 

2.75 

1.201 

-0.39044 

2.00 

2.0 

-0.41642 

3.00 

1.601 

-0.43342 

2.00 

2.2 

-0.41173 

3.00 

1.201 

-0.42923 

2.25 

1.6 

-0.44747 

3.25 

1.601 

-0.45831 

2.25 

1.8 

-0.44575 

3.25 

1.201 

-0.45554 

2.50 

1.2 

-0.44861 

3.50 

1.601 

-0.47547 

2.50 

1.6 

-0.47019 

3.50 

1.201 

-0.47310 

2.50 

1.8 

-0.46416 

4.00 

1.601 

-0.49496 

3.00 

1.2 

-0.48445 

4.00 

1.201 

-0.49229 

3.00 

1.6 

-0.49532 




geometry 

z zx 

geometry 

x — zx 

2.5 

1.8 

-0.40498 

2.5 

1.6 

-0.35863 

3.0 

1.2 

-0.45385 

2.5 

1.2 

—0.34837 

3.0 

1.6 

-0.46488 

3.0 

1.2 

-0.43792 

3.5 

1.2 

-0.48448 

3.0 

1.6 

-0.44049 

3.5 

1.6 

-0.49079 

3.5 

1.2 

-0.47799 




3.5 

1.6 

-0.47989 

geometry 

x - zy 

geometry 

X -f* zx 

2.5 

1.2 

-0.35527 

2.5 

1.6 

— 0.22158 

2.5 

1.6 

-0.37009 

2.5 

1.2 

— 0.21865 

3.0 

1.6 

-0.44546 

3.0 

1.6 

-0.37847 

3.0 

1.2 

-0.44091 

3.0 

1.2 

-0.37893 

3.5 

1.2 

-0.47941 

3.5 

1.2 

-0.45227 

3.5 

1.6 

-0.48221 

3.5 

1.6 

-0.45243 

geometry 

x, "h zy 


geometry 

x d - yx 

2.5 

1.6 

— 0.29103 

3.0 

1.6 

— 0.42983 
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2.5 

1.8 

-0.29236 

2.5 

2.0 

-0.29171 

2.5 

2.2 

-0.29083 

3.0 

1.2 

-0.39687 

3.0 

1.6 

-0.40586 

3.5 

1.2 

-0.45984 

3.5 

1.6 

-0.46419 

geometry 


z + yx 

2.2 

1.601 

-0.43573 

2.2 

1.801 

-0.43687 

2.2 

2.201 

-0.43270 

2.4 

1.601 

-0.45703 

2.4 

1.801 

-0.45402 

2.6 

1.601 

-0.47265 

2.6 

1.801 

-0.46658 

2.8 

1.201 

-0.46814 


3.0 

1.2 

-0.42693 

2.5 

1.2 

-0.32995 

2.5 

1.6 

-0.34194 

2.5 

1.8 

-0.33464 

3.5 

1.6 

—0.47373 

3.5 

1.2 

-0.47190 


geometry 


z + yx 

2.8 

1.601 

-0.48406 

3.0 

1.201 

-0.48000 

3.0 

1.601 

-0.49231 

3.2 

1.201 

-0.48839 

3.2 

1.601 

-0.49820 

3.5 

1.201 

-0.49648 

3.5 

1.601 

-0.50394 

4.0 

1.201 

-0.50299 

4.0 

1.601 

-0.50855 


a Distance between 0 and nearest H if tbe Hj bond points to the 0, otherwise the 
distance between 0 and the H 2 center of mass. 


27 


Table VI. Energies for H 2 4-H 2 0 4- 77 Eh with H 2 at its equilibrium geometry and 
H 2 0 displaced from equilibrium. All energies and distances are in atomic units. 


^ox 

gACPF 


-ft ox 

gACPF 



geometry 

z 4- zx, 

AS 


geometry 

z — ZX, 

AS 

3. 


-0.37285 


3. 


-0.43794 


4. 


-0.48282 


4. 


-0.49483 



geometry 

z + zx, 

SS 


geometry 

z — zx, 

SS 

3. 


— 0.36676 


3. 


-0.44101 


4. 


— 0.48045 


4. 


-0.49566 



geometry 

z + zx, 

B 


geometry 

z — zx, 

B 

4. 


-0.48928 


3. 


-0.44223 


3. 


-0.38958 


4. 


-0.49963 



a Distance between 0 and nearest H. 
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Table VII. Par am eters for non-bonding pair potentials in atomic units. 


A 

&ha 

cH A 

H 

1.618 

2.085 

0 

1.113 

1.384 
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Table VIII. Coefficients for bond-length dependence of H 2 long-range potential pa- 
rameters in atomic units. - 


goc 9u 92c n c 

- 1.26998 26.9001 9 0.882033 

- 0.176017 3.69209 9 0.715359 

og bj dq 

1.17467 1.72657 0.449685 7.26091 


r c 


2.85570 

1.8981 
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Table IX. Parameters for the vibrational dependence of the H 2 0 electrostatic prop- 


erties. See Eq. (33). 


M 

M° 00 

M 100 

M 010 

M no 

M 001 

Hz 

-1.034 

-0.1463 

-0.1463 

0.05971 

0.3387 

H* 

0.000 

0.1194 

-0.1194 

0.000 

0.000 

0° 

6.246 

-5.581 

-5.581 

3.307 

2.741 

e;. 

2.085 

-1.151 

-1.151 

0.2764 

-0.3180 

es. 

0.000 

1.552 

-1.552 

0.000 

0.000 

a 

-11.97 

7.545 

7.545 

-2.033 

0.3624 


° calculated with the origin at the center of mass of the H 2 0. 
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Table X. Angular expansion functions. 


no. 

P 

9 

m 

77l' 

-Apgmm' 

i 

0 

0 

0 

0 

l 

2 

0 

1 

0 

0 

cos /? 

3 

0 

2 

0 

0 

3 cos 2 /? — 1 

4 

2 

0 

0 

0 

3 cos 2 # — 1 

5 

2 

1 

0 

0 

cos/?( 3 cos 2 # — 1) 

6 

2 

2 

0 

0 

(3 cos 2 /? — l)( 3 cos 2 # — 1 ) 

7 

0 

2 

0 

2 

sin 2 /? cos(27) 

8 

2 

2 

2 

2 

sin 2 # cos 4 (/?/ 2 ) cos 2 (a 4 - 7) 

9 

2 

2 

2 

-2 

sin 2 # sin 4 (/?/2) cos 2(a — 7) 
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Table XI. Parameters for the fit of Eq.(43) in atomic units. 


p 

9 

m 

to' 

„sv 

v pgmm'l 

sv 

v SV 

u pqmm‘ 3 

0 

0 

0 

0 

8.0288800E— 01 

4.9208551E— 02 

— 1.9085039E— 03 

0 

1 

0 

0 

3.1244004E— 01 

— 9.0318002E— 03 

— 3.9700638E-04 

0 

2 

0 

0 

— 3.2853972E— 02 

1.7373376E— 03 

— 4.6648941E— 04 

2 

0 

0 

0 

— 1.2802147E— 01 

— 7.8669031E— 03 

-1.4491438E— 03 

2 

1 

0 

0 

-1.0989261E— 02 

1.8603093E— 03 

1.5477987E— 05 

2 

2 

0 

0 

1.3347056E— 02 

— 2.8587331E— 04 

— 6.2784176E— 05 

0 

2 

0 

2 

— 1.7133020E— 02 

5.8373975E— 04 

2.0297984E— 04 

2 

2 

2 

2 

6.0174056E— 02 

2.5171433E— 04 

— 6.3496729E— 05 

2 

2 

2 

-2 

— 1.4148450E— 01 

1.5470255E— 02 

— 3.9044009E— 04 

r 5fl 

d SR 

d 

€ 

a* 

r* 


0. 

1.72 

4.0 

0.562 

2.5 

3.0 
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Table XII. Parameters for the fit of Eq.(44) in atomic units. 


p 


m 

m' 

r.sv 

v SV 

pqmm 1 2j 

- sv 

Vpqmm 1 3 j 

0 

0 

0 

0 

j = 1 

— 7.8409154E— 01 

5.7772398E— 02 

— 5.2699265E— 03 

0 

1 

0 

0 

— 2.8981900E— 02 

— 6 .4440092E — 03 

1.8197968E— 04 

0 

2 

0 

0 

— 2.7464954E— 02 

5.4274961E— 03 

7.8214372E— 05 

2 

0 

0 

0 

— 1.7779328E— 01 

— 2.2780539E— 02 

— 2.5086883E— 04 

2 

1 

0 

0 

— 3.3807402E— 02 

1.3647800E— 03 

— 4.6390793E— 04 

2 

2 

0 

0 

2.2092459E— 02 

-8.8503368E— 04 

4.5817703E— 04 

0 

2 

0 

2 

2.6905372E— 02 

— 8 .4002092E— 03 

2.6707433E— 03 

2 

2 

2 

2 

8.2026964E— 02 

1.2664571E— 04 

2.8665551E— 05 

2 

2 

2 

-2 

— 2.3292928E— 01 
j = 2 

2.3490078E— 02 

— 2.2725534E— 03 

0 

0 

0 

0 

— 2.6892232E— 01 

— 3.7369693E— 02 

2.9057342E— 03 

0 

1 

0 

0 

1.7217490E— 01 

— 1.5812969E— 02 

1.5763320E— 03 

0 

2 

0 

0 

1.1942930E— 01 

— 1 .2523267E— 02 

2.4867362E— 03 

2 

0 

0 

0 

— 2.0629499E— 01 

— 1.1532656E— 02 

— 5.2182780E— 04 

2 

1 

0 

0 

1.8435082E— 01 

— 2.7919859E— 02 

4.9524319E— 03 

2 

2 

0 

0 

1.0729963E— 01 

— 6.8764652E— 03 

1.2178887E— 03 

0 

2 

0 

2 

— 4.1360738E— 02 

3.7667774E— 03 

— 7.1684775E— 04 

2 

2 

2 

2 

1.4759330E— 02 

1.0189559E— 03 

— 9.5442150E— 05 

2 

2 

2 

-2 

— 3.7665040E— 02 

— 2.4012130E— 03 

1.2519036E— 03 
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Table XIII. Comparison of analytical representation to the ab initio data of Table 


VI. All energies are in mEh- 


r ( 00 ) 

H 2 0 6 

geometry 

yint 

ab initio 

fit 

3. 

B 

x + zx 

122.1 

114.1 

3. 

ss 

X + zx 

139.9 

139.2 

3. 

AS 

X + zx 

133.6 

134.4 

4. 

B 

X + zx 

22.41 

21.08 

4. 

SS 

X + zx 

26.25 

24.36 

4. 

AS 

X + zx 

23.66 

23.61 

3. 

B 

x — ZX 

69.45 

70.10 

3. 

SS 

x — ZX 

65.69 

67.25 

3. 

AS 

x — zx 

68.54 

68.88 

4. 

B 

x — zx 

12.05 

11.99 

4. 

SS 

x — zx 

11.04 

11.20 

4. 

AS 

x — zx 

11.66 

11.61 


k H 2 0 geometry. 
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Figure captions: 

Fig. 1: Comparison of V 0 to the ob initio data of Table III. Both have 0.01 Eh 
added to them. The solid line and o are for the z + zx orientation and the dashed 
line and a are for the z — zx configuration. 

Fig. 2: Comparison of the fit to the H 2 polarizability to the ab initio data from 
Ref. 25. The solid line and o are for the || component and dashed line and □ are for 
the -L component. 

Fig. 3: Comparison of V’ nt to data for z + yx orientation. The variable r changes 
with r_ to maintain a fixed Roh distance. The curves are for J2oh=2.2, 2.4, 2.6, 
2.8, 3.0, 3.2, 3.5, and 4.0 a 0 . 
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Theoretical Characterization of the Potential Energy Surface 

for 

H + O 2 —> ho; -> HO + 0 

III. Computed Points to Define a Global Potential for H + O2 

Stephen P. Walch a and Ronald J. Duchovic 0 
ELORET Institute 
Sunnyvale, Ca. 94087 

Abstract. Recent calculations on the H + 0 2 surface have focused on the minimum 
energy path region of the surface ( J. Chem. Phys., 88, 6273(1988) I ) and on 
the saddle point region for H atom exchange via a T-shaped H0 2 complex ( J. 
Chem. Phys., 91, 2373(1989) II ). In this paper, additional computed points are 
reported which, when combined with previously reported points, permit a global 
representation of the H + 0 2 and H0 2 regions of the potential energy surface. The 
calculations axe complete active space SCF/ externally contracted configuration 
interaction (CASSCF/CCI) using the same wavefunction as in II. The new points 
characterize the potential for all angles of approach ranging from perpendicular to 
collinear with the 00 bond, for H to center of mass of 0 2 distances ranging from 
5.0 ao in to a distance corresponding to greater than 30 kcal/mol up the inner wall. 
A new collinear exchange saddle point is reported. 

“Mailing Address: NASA Ames Research Center, Moffett Field, CA 94035. 
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I. Introduction 


The reaction 


H + 0 2 -* HO; -* HO + 0 (1) 

is an important reaction in combustion. In addition, as part of a program to model 
combustion processes important in the design of the high-speed civil transport, rate 
constants for H atom diffusion in air are needed. Computation of these rate con- 
stants requires a global potential for H + O 2 . Recently, calculated points have been 
reported for the minimum energy path (MEP) region of the potential energy surface 
for reaction (1) [1]. (Hereafter, referred to as I.) Also, additional calculations have 
been reported which characterize the pathway for exchange of an H atom of HO 2 
via a T-shaped HO 2 saddle point [2]. (Hereafter, referred to as II.) In the present 
paper, additional points are computed which, together with the points previously 
computed in II, provide the data needed to produce a global representation of the 
H + O 2 interaction. 

The computational method is discussed in Section II, the results are presented in 
Section III, and the conclusions are given in Section IV. 

II. Computational Details. 

The coordinate system for these calculations is the same as that used in II. The 
coordinates used are the 00 distance (r oo)> the H to center of mass of O 2 distance 
(rn-oo)j and the angle (6) between a line connecting H to the center of mass of O 2 
and a normal to the 00 bond at the bond midpoint ( 6 is 0° for T-shaped H-O 2 
and 90° for collinear H-02). 

The calculations were carried out in the same way as in I. A problem which was 
encountered in designing the calculations in II was that in the vicinity of the T- 
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shaped HO 2 saddle point a 5a'2a" CASSCF active space was needed while in the 
remainder of the surface a 4a'la" active space was more appropriate. The way in 
which this was handled was to use the larger active space for 0 = 0°, 10°, and 20° 
and the smaller active space for larger 0 values. In the current calculations the 
smaller active space was used for all points except those having 0 = 0°, 10°, and 
20° and th-oo smaller than the values considered in II. The orbitals obtained using 
these two active space sizes are similar enough that there has been no indication 
of discontinuities in the subsequent CCI energies which are CASSCF/CCI with a 
5a' 2a" CASSCF wavefunction as reference space for the Cl. All of the calculations 
were carried out in C, symmetry. Some of the calculations in II were carried out in 
Cj 2 v sy mm etry and where these were combined with the present calculations this is 
indicated in the tables. 

Most of the calculations were carried out with the [4s3p2dlf/3s*2pld] basis set 
described in II. As discussed in II, the notation s* indicates that the contraction of 
the first natural orbital for H is based on the atomic SCF orbital, while in the original 
basis of Almloff and Taylor [3] this contraction is based on natural orbitals from 
a Cl calculation on H 2 . Additional calculations were carried out along the MEP 
defined with the [4s3p2dlf/3s*2pld] basis set using a larger [5s4p3d2f/4s*3pd2d] 
basis set [3]. 

Calculations were carried out for 0 values of 0° through 90° in 10° increments. 
For each 0 value a grid was carried out over r h-oo and T oo to obtain minimum 
energy cuts at fixed 0 values. 

The calculations were carried out on the NASA Ames Cray Y-MP/832. These 
calculations used the MOLECULE[4]-SWEDEN[5] system of programs. 

III. Results and Discussion 

The computed energies are given in Table AI of the appendix. In order to aid 
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in visualizing the surface, the energy was evaluated along fixed 9 minimum energy 
cuts. For each 9 and r h-oo, roo was varied and the energy at the minimum 
and the optimal too are given in Table I. These minimum energy cuts are shown 
graphically in Figs. 1 and 2. 

A number of features of the surface are evident in Figs. 1 and 2. These features 
include a very small entrance channel barrier, a minimum ~ 51 kcal/mol below H + 
O 2 , a T-shaped(0 = 0 °) H exchange saddle point « 13 kcal/mol below H + O 2 , and 
a collinear(0 = 0 °) H exchange saddle point ~ 9 kcal/mol above H 4- O 2 . It is also 
evident from Figs. 1 and 2 that for each 9 value there is a maximum and subsequent 
minimum as th-oo is decreased. This leads to additional stationary points which 
are maxima on the global surface for 9 = 0° and 90°. In addition to these features 
there is also am OH + 0 product channel which is discussed below. The MEP for 
this process corresponds approximately to increasing too while varying the other 
two geometrical parameters. 

In order to define the MEP for the reactants channel, polynomial fits ( six-term 
quadratic in too and 9 ) were obtained at each th-oo distance using the three 
9 values nearest the minimum and three values of too (9 points). At the HO 2 
minimum a ten-term quadratic polynomial in sill three coordinates was obtained 
to define the stationary points. These results are given in Table II. Additional 
calculations were carried out along this MEP using the [5s4p3d2f/4s*3p2d] basis 
set. These results are also given in Table II. 

In Fig. 3, the energy along the MEP with the [4s3p2dlf/3s*2pld] basis set is 
compared with the energy along the MEP of I (using the basis set of I). Both 
curves are referenced to the respective H + O 2 asymptotic energy, therefore they 
coincide at large th-oo • From Fig. 3 it is seen that the curves essentially coincide 
up to th-oo « 3.5a<). For shorter th-oo the comparison is more difficult since 
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the MEP of I was determined on a courser th-oo grid than in the present work. 
For I only the energy at the minimum is shown near the bottom of the well. The 
saddle point geometries and relative energies for the H-O 2 entrance channel saddle 
point and HOj minimum are given in Table III. Here it is seen that the geometry 
at the HO 2 minimum is in good agreement with that obtained in I (maximum 
error is 0.02 ao in toh ) hut the H-O 2 well depth is 1.9 kcal/mol smaller with the 
[4s3p2dlf/3s*2pld] basis set and 0.2 kcal/mol smaller with the [5s4p3d2f/4s*3p2d] 
basis set. The H-O 2 entrance saddle point geometry differs more between the two 
calculations (maximum error 0.18 ao in r/jo, but the difference in the barrier height 
is only 0.04 kcal/mol for the larger basis set). 

Using the computed HO 2 well depth, the computed HO 2 and O 2 vibrational 
frequencies from I, and the experimental H atom AH® of 51.6 kcal/mol [6], gives 
a AH® of 5.3 kcal/mol for HO 2 . This value may be compared to the value of 
3.5^5 kcal/mol recommended by Benson [7]. Using the error limits quoted by 
Benson, the error in the computed HO 2 well depth is between 0.8 and 2.3 kcal/mol. 
This error range is quite reasonable for the basis set and level of calculation used 
here, and thus these calculations are consistent with the AH® for HO 2 recommended 
by Benson. 

Fig. 4 and Fig. 5 compare the variation of 6 and r 00 along the MEP as a function 
of r h-oo for the present work and the calculations in I. The general features of the 6 
and too variations are similar but there are some significant variations in detail. It 
is interesting that these geometrical parameters seem to coincide at the stationary 
points but vary in between. This may result from the non-rigorous method of 
defining the MEP which is used here. It should be noted that the MEP as defined 
here is mainly for purposes of visualizing the surface. A more accurate definition of 
these features must await the development of an accurate analytic representation 
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of the surface. 

Table IV shows three fixed 9 minimum energy cuts in which rjf-oo is varied 
for each too • These cuts are for motion orthogonal to the reactant channel and 
represent the intersection of the product channel MEP and the H00 minimum. 
These curves are plotted in the inset to Fig. 1. 

Fig. 1 also shows higher energy regions of the surface corresponding to 9 values 
near 0° and 90°. It is evident from Fig. 1 that there are saddle points for H atom 
exchange along cuts with 0 = 0° and 90°. For these two choices of 6 the actual 
symmetry is C 20 and the barrier (maximum on the surface) prior to the saddle 
point corresponds to a curve crossing ( 2 i?i —* 2 .4.2 for 9 = 0° and 2 E - — * 2 II for 9 
= 90°). Table III also shows the saddle point geometries and barriers for the two 
exchange saddle points. Here it is seen that the T-shaped exchange saddle point is 
13 kcal/mol below H + O 2 , while the collinear saddle point is 9 kcal/mole above 
H + O 2 . Both of these barriers are below the OH + O product channel. Given 
the barriers to dissociation for 9 = 90° and also possible centrifugal barriers, it is 
possible that a complete circular motion of H around O 2 could occur for energies 
greater than 9 kcal/mol. In any case, the T-shaped H exchange saddle point is 
accessible to all H + O 2 collisions, since it is below the H + O 2 asymptote. It is 
anticipated that these previously uncharacterized features of the H + O 2 potential 
will have significant dynamical consequences. 

IV. Conclusions. 

Computed points are reported which, when combined with previously published 
points from II, permit a global representation of the H + O 2 and HO 2 regions of 
the potential energy surface for the reaction: 

h + o 2 -► ho; -» HO + 0 
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The new points characterize the potential for all angles of approach of H to O 2 
and for the inner repulsive wall region. In addition to connecting the points for 
the T-shaped H-O 2 exchange saddle point ( previously characterized in II ) to the 
reactant minimum energy path and inner wall regions, this work also characterizes 
a new collinear H-O 2 exchange saddle point which is only 9 kcal/mol above H + 
O 2 . It is anticipated that these new features of the potential energy surface will 
have significant dynamical consequences. 
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Table I . H + O2 Energy Along Fixed 0 Cuts . 


TH-00 

O 

1 

O 

e 

Energy 

5E 

1.4 

2.878 

0.0 

- 0.46948 

- 3.56 

1.2 

3.101 

0.0 

- 0.44941 

9.04 

5.0 

2.28 

10.0 

- 0.46208 

1.09 

4.5 

2.28 

10.0 

- 0.45963 

2.62 

4.0 

2.28 

10.0 

- 0.45435 

5.94 

3.5 

2.283 

10.0 

- 0.44366 

12.64 

3.0 

2.335 

10.0 

- 0.42631 

23.53 

2.5 

2.603 

10.0 

- 0.45009 

8.61 

2.0 

2.627 

10.0 

- 0.48636 

- 14.15 

1.8 

2.670 

10.0 

- 0.49069 

- 16.87 

1.6 

2.745 

10.0 

- 0.48507 

- 13.34 

1.4 

2.903 

10.0 

-.46851 

- 2.85 

1.2 

3.164 

10.0 

-.44622 

11.04 

5.0 

2.28 

20.0 

- 0.46254 

0.80 

4.5 

2.28 

20.0 

- 0.46078 

1.90 

4.0 

2.28 

20.0 

- 0.45737 

4.04 

3.5 

2.300 

20.0 

- 0.45235 

7.19 

3.0 

2.412 

20.0 

- 0.45429 

5.97 

2.5 

2.545 

20.0 

- 0.48478 

- 13.16 

2.2 

2.509 

20.0 

- 0.50633 

- 26.68 
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2.0 

2.594 

20.0 

- 0.50976 

- 28.83 

1.8 

2.664 

20.0 

- 0.50461 

- 25.60 

1.6 

2.802 

20.0 

- 0.48671 

- 14.37 

1.4 

3.060 

20.0 

- 0.45957 

2.66 

1.2 

3.653 

20.0 

- 0.43836 

15.97 


5.0 

2.28 

30.0 

- 0.46308 

0.46 

4.5 

2.28 

30.0 

- 0.46227 

0.97 

4.0 

2.284 

30.0 

- 0.46168 

1.34 

3.5 

2.324 

30.0 

- 0.46479 

- 0.61 

3.0 

2.444 

30.0 

- 0.48378 

- 12.5 

2.8 

2.494 

30.0 

- 0.49779 

- 21.3 

2.6 

2.521 

30.0 

- 0.51274 

- 30.7 

2.4 

2.543 

30.0 

- 0.52515 

- 38.5 

2.2 

2.555 

30.0 

- 0.53097 

- 42.1 

2.0 

2.557 

30.0 

- 0.52470 

- 38.2 

1.8 

2.666 

30.0 

- 0.50064 

- 23.1 

1.6 

2.974 

30.0 

- 0.45961 

2.64 


5.0 

2.28 

40.0 

- 0.46329 

0.33 

4.5 

2.28 

40.0 

- 0.46315 

0.41 

4.0 

2.294 

40.0 

- 0.46499 

- 0.74 

3.5 

2.36 

40.0 

- 0.47546 

- 7.31 

3.0 

2.487 

40.0 

- 0.50753 

- 27.4 

2.8 

2.516 

40.0 

- 0.52402 

- 37.8 
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2.6 

2.530 

40.0 

- 0.53803 

- 46.6 

2.4 

2.528 

40.0 

- 0.54451 

- 50.6 

2.2 

2.523 

40.0 

- 0.53653 

- 45.6 

2.0 

2.533 

40.0 

- 0.50476 

- 25.7 

1.8 

2.648 

40.0 

- 0.43845 

15.9 

5.0 

2.28 

50.0 

- 0.46262 

0.75 

4.5 

2.28 

50.0 

- 0.46208 

1.09 

4.0 

2.357 

50.0 

- 0.46450 

- 0.43 

3.5 

2.416 

50.0 

- 0.48006 

- 10.2 

3.0 

2.435 

50.0 

- 0.52048 

- 35.6 

2.8 

2.537 

50.0 

- 0.53564 

- 45.1 

2.6 

2.515 

50.0 

- 0.54252 

- 49.4 

2.4 

2.472 

50.0 

- 0.53321 

- 43.5 

2.2 

2.405 

50.0 

- 0.49595 

- 20.2 

2.0 

2.303 

50.0 

- 0.41303 

31.9 

5.0 

2.28 

60.0 

- 0.46076 

1.91 

4.5 

2.28 

60.0 

- 0.45806 

3.61 

4.0 

2.313 

60.0 

- 0.45677 

4.42 

3.5 

2.491 

60.0 

- 0.47462 

- 6.78 

3.0 

2.559 

60.0 

- 0.51669 

- 33.2 

2.5 

2.404 

60.0 

- 0.50224 

- 24.1 

2.2 

2.252 

60.0 

- 0.39683 

42.0 

5.0 

2.28 

70.0 

- 0.45805 

3.61 


11 



4.5 

2.28 

70.0 

- 0.45152 

7.71 

4.0 

2.297 

70.0 

- 0.44194 

13.72 

3.5 

2.588 

70.0 

- 0.45769 

3.84 

3.0 

2.543 

70.0 

- 0.49354 

- 18.66 

2.5 

2.296 

70.0 

- 0.43458 

18.34 


5.0 

2.28 

80.0 

- 0.45545 

5.25 

4.5 

2.28 

80.0 

- 0.44514 

11.72 

4.0 

2.255 

80.0 

- 0.42539 

24.11 

3.5 

2.674 

80.0 

- 0.43334 

19.12 

3.0 

2.503 

80.0 

- 0.46114 

1.68 

2.5 

2.222 

80.0 

- 0.37282 

57.1 


5.0 

2.28 

90.0 

- 0.45456 

5.80 

4.5 

2.28 

90.0 

- 0.44269 

13.25 

4.0 

2.28 

90.0 

- 0.41797 

28.77 

3.5 

2.723 

90.0 

- 0.41288 

31.96 

3.25 

2.617 

90.0 

- 0.43737 

16.6 

3.0 

2.493 

90.0 

- 0.44634 

10.96 

2.75 

2.352 

90.0 

- 0.42450 

24.7 
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Table II. Computed MEP for H Atom Addition to O 2 . 


r H—00 

ro-o 

6 

Energy(bs 1°) 

Energy(bs 2 b ) 


10.0 

2.28 


(-0.46381) 

0.0 -150.46219(-.47133) 

0.0 

5.0 

2.28 

37.4 

(-0.46332) 

0.31 -150.46134(-.47068) 

0.41 

4.60 

2.28 

40.6 

(-0.46325) 

0.35 -150.46110(-.47060) 

0.46 

4.5 

2.28 

39.5 

(-0.46315) 

0.41 


4.0 

2.30 

43.2 

(-0.46520) 

-0.87 -150.46287(-.47304) 

-1.07 

3.5 

2.41 

49.6 

(-0.48008) 

-10.21 -150.47626(-.48805) 

-10.49 

3.0 

2.54 

52.7 

(-0.52111) 

-35.96 -150.51830(-.53019) 

-36.9 

2.8 

2.55 

52.8 

(-0.53612) 

-45.38 -150.53274(-. 54445) 

-45.9 

2.6 

2.52 

47.2 

(-0.54331) 

-49.89 -150.54244(-. 55393) 

-51.8 

2.49 

2.52 

43.5 

(-0.54546) 

-51.24 -150.54415(-. 55557) 

-52.9 

2.4 

2.51 

41.3 

(-0.54462) 

-50.71 -150.54393(-. 55529) 

-52.7 

2.2 

2.52 

36.1 

(-0.53994) 

-47.77 -150.53760(-.54886) 

-48.7 


0 [4s3p2dlf/3s*2p] ANO basis set. 

6 [5s4p3d2f/4s*3p2d] ANO basis set. 
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Table III. Stationary Points on the H + O 2 Surface. 



H-0 2 

O 

M 

H-O 2 (0=O°) c 

H-O2(0: 

TOO 

2.28 

2.52 

2.70 

2.51 

r H-00 

4.60 

2.49 

1.74 

3.06 

e 

46.3 

43.5 

0.0 

90.0 

toh 

3.96 

1.86 



gla 

118 

104.2 



E(bs 1) 

0.35 

-51.2 

-13 

9.2 

E(bs 2) 

0.46 

-52.9 



T b 

T oo 

2.29 

2.52 



T b 

t oh 

4.14 

1.84 



Q,b 

116.4 

104.4 



E fr 

0.5 

-53.1 




a L HOO. 
b values from ref. I. 
c from ref. II. 
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Table IV. H + O 2 Energy Along Fixed 0 Cuts. 


T'o-o 

TH-OO 

e 

Energy 

£E 

2.4 

2.410 

40.0 

-0.54100 

-48.4 

2.5 

2.410 

40.0 

-0.54436 

-50.4 

2.6 

2.412 

40.0 

-0.54341 

-49.9 

2.7 

2.415 

40.0 

-0.53944 

-47.5 


2.4 

2.228 

30.0 

-0.52560 

-38.8 

2.5 

2.212 

30.0 

-0.53031 

-41.7 

2.6 

2.199 

30.0 

-0.53050 

-41.8 

2.7 

2.176 

30.0 

-0.52758 

-40.0 


2.4 

2.599 

50.0 

-0.53967 

-47.6 

2.5 

2.614 

50.0 

-0.54251 

-49.4 

2.6 

2.630 

50.0 

-0.54116 

-48.5 

2.7 

2.646 

50.0 

-0.53693 

-45.9 
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Figure Captions. 


Fig. 1. Potential surface for H + O 2 . The figure shows ten fixed 6 minimum energy 
cuts. For each r h-oo> roo was varied and the minimum energy is shown in the 
figure. These cuts are for the reactant (H + 0 2 ) channel. The inset shows three 
fixed 6 minimum energy cuts. Here r h—oo is varied for each too- These cuts are 
for motion orthogonal to the reactant channel and represent the intersection of the 
product channel MEP and the H00 minimum. 

Fig. 2. Potential surface for H + 0 2 . This figure shows the same information as 
Fig. 1 in the format of a perspective plot. 

Fig. 3. Comparison of energy as a function of tjj—oo along the MEP from I and 
the present calculations. For the present calculations results are shown with both 
basis sets. 

Fig. 4. Comparison of 6 as a function of th—oo for the MEP from I and the present 
work. 

Fig. 5. Comparison of too as a function of rjj— 00 f° r the MEP from I and the 
present work. 
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THE H + 0« SURFACE 






Walch 




"H0©"* 


□ [4s3p2d1f/3s*2p] 

O [5s4p3d2f/4s*3p2d] 






Appendix. The appendix contains a table of all the computed CASSCF/CCI ener- 
gies. The energies are in the form CCI(CCI +Q). Note that for the CCI + Q energies 
-150. is not repeated. Thus, for the first point the CCI energy is -150.45491 and 
the CCI+Q energy is -150.46381. Distances are in ao, angles are in degrees, and 
energies axe in Eh. 
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Table AI. H + O 2 Long Range Points. 


TH-OO 

O 

1 

O 

e 

10.0 

2.28 


11.15 

2.3 

90.0 

6.15 

2.3 

90.0 

5.15 

2.3 

90.0 

4.15 

2.3 

90.0 

5.0 

2.28 

90.0 

4.5 

2.28 

90.0 

4.0 

2.28 

90.0 

3.5 

2.4 

90.0 

3.5 

2.5 

90.0 

3.5 

2.6 

90.0 

3.5 

2.7 

90.0 

3.5 

2.8 

90.0 

3.25 

2.5 

90.0 

3.25 

2.6 

90.0 

3.25 

2.7 

90.0 

3.0 

2.4 

90.0 

3.0 

2.5 

90.0 

3.0 

2.6 

90.0 

3.0 

2.7 

90.0 

2.75 

2.3 

90.0 

2.75 

2.4 

90.0 


Energy 


- 150 . 45491 (-. 46381 ) 

- 150 . 45482 (-. 46385 )° 
- 150 . 45364 (-. 46269 )° 
- 150 . 44734 (- . 45656 )“ 
- 150 . 41740 (-. 42687 )° 

- 150 . 44542 (-. 45456 ) 
- 150 . 43341 (-. 44269 ) 
- 150 . 40842 (- . 41797 ) 

- 150 . 38260 (-. 39467 ) 
- 150 . 39275 (-. 40517 ) 
- 150 . 39803 (-. 41083 ) 
- 150 . 39960 (-. 41281 ) 
- 150 . 39837 (-. 41206 ) 

- 150 . 42224 (-. 43441 ) 
- 150 . 42481 (-. 43731 ) 
- 150 . 42299 (-. 43588 ) 

- 150 . 43172 (-. 44341 ) 
- 150 . 43434 (-. 44632 ) 
- 150 . 43021 (-. 44251 ) 
- 150 . 42044 (-. 43310 ) 
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- 150 . 41175 (- . 42306 ) 
- 150 . 41165 (- . 42325 ) 



2.75 

2.5 

90.0 

-150.40087(-.41273) 

2.75 

2.6 

90.0 

-150.38048(-. 39260) 

5.0 

2.28 

80.0 

-150.44640(-. 45545) 

4.5 

2.28 

80.0 

-150.43592(-.44514) 

4.0 

2.28 

80.0 

-150.41553(-.42522) 

4.0 

2.3 

80.0 

-150.41504(- .42483) 

4.0 

2.4 

80.0 

-150.40886(-.41950) 

3.5 

2.28 

80.0 

-150.39797(-.41069) 

3.5 

2.3 

80.0 

-150.40051(-.41324) 

3.5 

2.4 

80.0 

-150.41073(-.42345) 

3.5 

2.5 

80.0 

-150.41696(-.42971) 

3.5 

2.6 

80.0 

-150.41991(-.43278) 

3.5 

2.7 

80.0 

-150.42030(-.43327) 

3.5 

2.8 

80.0 

-150.41872(-. 43168) 

3.0 

2.4 

80.0 

-150.44633(-.45797) 

3.0 

2.5 

80.0 

-150.44927(-.46114) 

3.0 

2.6 

80.0 

-150.44627(-.45833) 

3.0 

2.7 

80.0 

-150.43848(-.45071) 

2.5 

2.1 

80.0 

-150.34869(-.35939) 

2.5 

2.2 

80.0 

-150.36147(-.37237) 

2.5 

2.3 

80.0 

-150.35636(-. 36744) 

2.5 

2.4 

80.0 

-150.33626(-. 34751) 

2.5 

2.5 

80.0 

-150.30295(-.31437) 

5.0 

2.28 

70.0 

-150.44899(-. 45805) 

4.5 

2.28 

70.0 

-150.44221(-. 45152) 
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4.0 

2.28 

70.0 

-150.43188(-.44187) 

4.0 

2.3 

70.0 

-150.43181(-.44194) 

4.0 

2.4 

70.0 

-150.42845(- .43937) 

3.5 

2.28 

70.0 

-150.43284(-.44477) 

3.5 

2.3 

70.0 

-150.43473(-.44677) 

3.5 

2.4 

70.0 

-150.44153(- .45387) 

3.5 

2.5 

70.0 

-150.44479(-.45720) 

3.5 

2.6 

70.0 

-150.44536(-.45768) 

3.5 

2.7 

70.0 

-150.44456(-.45689) 

3.0 

2.4 

70.0 

-150.47718(-.48871 ) 

3.0 

2.5 

70.0 

-150.48143(-.49310) 

3.0 

2.6 

70.0 

-150.48094(-.49270) 

3.0 

2.7 

70.0 

-150.47690(-.48865) 

2.5 

2.2 

70.0 

-150.41862(-.4291 1 ) 

2.5 

2.3 

70.0 

-150.42388(-.43457) 

2.5 

2.4 

70.0 

-150.41739(-.42823) 

2.5 

2.5 

70.0 

-150.40149(-.41242) 

5.0 

2.28 

60.0 

-150.45167(-.46076) 

4.5 

2.28 

60.0 

-150.44871(-.45806) 

4.0 

2.28 

60.0 

-150.44649(-.45649) 

4.0 

2.3 

60.0 

-150.44661(-.45673) 

4.0 

2.4 

60.0 

-150.44399(-.45473) 

3.5 

2.28 

60.0 

-150.45585(-.46717) 

3.5 

2.3 

60.0 

-150.45732(-.46875) 

3.5 

2.4 

60.0 

-150.46175(-.47358) 

3.5 

2.5 

60.0 

-150.46258(-.47461) 
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3.5 

2.6 

60.0 

-150.46107(-.47310) 

3.0 

2.3 

60.0 

-150.49059(-. 50188) 

3.0 

2.4 

60.0 

-150.50007(-. 51146) 

3.0 

2.5 

60.0 

-150.50446(-.51597) 

3.0 

2.6 

60.0 

-150.50480(-. 51634) 

2.5 

2.3 

60.0 

-150.48774(-.49838) 

2.5 

2.4 

60.0 

-150.49144(-. 50223) 

2.5 

2.5 

60.0 

-150.48810(-.49898) 

2.2 

2.2 

60.0 

-150.38533(-. 39542) 

2.2 

2.3 

60.0 

-150.38534(-.39559) 

2.2 

2.4 

60.0 

-150.37485(-.38521 ) 

2.2 

2.5 

60.0 

-150.35665(-.36707) 

5.0 

2.28 

50.0 

-150.45354(-.46262) 

4.5 

2.28 

50.0 

-150.45276(-.46208) 

4.0 

2.28 

50.0 

-150.45418(-.46407) 

4.0 

2.3 

50.0 

-150.45427(-. 46425) 

4.0 

2.4 

50.0 

-150.45131(-. 46179) 

4.0 

2.5 

50.0 

-150.44439(- .45534) 

4.0 

2.6 

50.0 

-150.43532(-.44663) 

3.5 

2.28 

50.0 

-150.46541(-.47630) 

3.5 

2.3 

50.0 

-150.46646(-.47745) 

3.5 

2.4 

50.0 

-150.46860(- .48001 ) 

3.5 

2.5 

50.0 

-150.46703(-.47869) 

3.5 

2.6 

50.0 

-150.46319(-.47493) 

3.0 

2.3 

50.0 

-150.49815(*. 50933) 
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3.0 

2.4 

50.0 

-150.50602(-.51735) 

3.0 

2.5 

50.0 

-150.50888(-.52027) 

3.0 

2.6 

50.0 

-150.50834(-.51976) 

2.8 

2.4 

50.0 

-150.52075(-.53186) c 

2.8 

2.5 

50.0 

-150.52416(-.53537) c 

2.8 

2.6 

50.0 

-150.52359(-.53484) e 

2.6 

2.4 

50.0 

-150.52875(-.53967) c 

2.6 

2.5 

50.0 

-150.53146(-.54247) c 

2.6 

2.6 

50.0 

-150.52991(-.54096) c 

2.4 

2.4 

50.0 

-150.52137(-.53206) c 

2.4 

2.5 

50.0 

-150.52225(-.53304) c 

2.4 

2.6 

50.0 

-150.51875(-.52958) c 

2.2 

2.3 

50.0 

-150.48128(-.49161) 

2.2 

2.4 

50.0 

-150.48548(-.49595) c 

2.2 

2.5 

50.0 

- 1 50.48344(- .49399 ) c 

2.2 

2.6 

50.0 

-150.47714(-.48772) c 

2.0 

2.3 

50.0 

-150.40079(-.41089) 

2.0 

2.4 

50.0 

-150.40111(-.41133) c 

2.0 

2.5 

50.0 

-150.39575(-.40604) c 

2.0 

2.6 

50.0 

-150.38684(-.39716) c 

1.8 

2.6 

50.0 

-150.22628(-.23634) c 

1.8 

2.5 

50.0 

-150.23443(-.24445) c 

1.8 

2.4 

50.0 

-150.24079(-.25074) c 

1.6 

2.6 

50.0 

-149.98075(-.99050) c 

1.6 

2.5 

50.0 

-149.97765(-.98735) c 
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1.6 

2.4 

50.0 

-149.97579(-.98544) c 

5.0 

2.28 

40.0 

-150.45423(-.46329) 

4.5 

2.28 

40.0 

-150.45391(-.46315) 

4.0 

2.28 

40.0 

-150.45524(-.46493) 

4.0 

2.3 

40.0 

-150.45521(-.46498) 

4.0 

2.4 

40.0 

-150.45152(-.46173) 

4.0 

2.5 

40.0 

-150.44363(-.45425) 

4.0 

2.6 

40.0 

-150.43332(-.44430) 

3.5 

2.28 

40.0 

-150.46344(-.47399) 

3.5 

2.3 

40.0 

-150.46409(-. 47473) 

3.5 

2.4 

40.0 

-150.46411(-.47515) 

3.5 

2.5 

40.0 

-150.46028(-.47159) 

3.5 

2.6 

40.0 

-150.45418(-.46562) 

3.0 

2.3 

40.0 

-150.48923(-. 50028) 

3.0 

2.4 

40.0 

-150.49505(-. 50626) 

3.0 

2.5 

40.0 

-150.49619(-. 50750) 

3.0 

2.6 

40.0 

-150.49407(-. 50541) 

2.8 

2.4 

40.0 

-150.51051(-.52162) c 

2.8 

2.5 

40.0 

-150.51279(-.52397) c 

2.8 

2.6 

40.0 

-150.51153(-.52276) c 

2.6 

2.4 

40.0 

-150.52366(-.53460) c 

2.6 

2.5 

40.0 

-150.52680(-.53785) c 

2.6 

2.6 

40.0 

-150.52592(-.53704) c 

2.4 

2.4 

40.0 

-150.53019(-.54098) c 

2.4 

2.5 

40.0 

-150.53343(-.54434) c 
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2.4 

2.6 

40.0 

-150.53241(-.54338) c 

2.2 

2.4 

40.0 

-150.52262(-.53324) c 

2.2 

2.5 

40.0 

-150.52569(-.53641) c 

2.2 

2.6 

40.0 

-150.52448(-.53526) c 

2.0 

2.4 

40.0 

-150.49070(-.50112) c 

2.0 

2.5 

40.0 

-150.49402(-.50454) c 

2.0 

2.6 

40.0 

-150.49326(-.50384) c 

1.8 

2.8 

40.0 

-150.42584(-.43624) 

1.8 

2.6 

40.0 

-150.42786(-.43823) c 

1.8 

2.5 

40.0 

-150.42605(-.43636) c 

1.8 

2.4 

40.0 

-150.42065(-.43086) c 

1.6 

2.6 

40.0 

-150.31825(-.32841) c 

1.6 

2.5 

40.0 

-150.30895(-.31903) c 

1.6 

2.4 

40.0 

-150.29676(-.30674) c 

5.0 

2.28 

30.0 

-150.45406(-.46308) 

4.5 

2.28 

30.0 

-150.45312(-.46227) 

4.0 

2.28 

30.0 

-150.45219(-.46167) 

4.0 

2.3 

30.0 

-150.45205(-.46160) 

4.0 

2.4 

30.0 

-150.44768(-.45760) 

4.0 

2.5 

30.0 

-150.43888(-.44920) 

4.0 

2.6 

30.0 

-150.42747(-.43815) 

3.5 

2.28 

30.0 

-150.45401(- .46423) 

3.5 

2.3 

30.0 

-150.45431(-.46462) 

3.5 

2.4 

30.0 

-150.45244(- .46317) 

3.5 

2.5 

30.0 

- 1 50 .44664(- .45766 ) 
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3.5 

2.6 

30.0 

-150.43856(-. 44975) 

3.0 

2.3 

30.0 

-150.46833(-.47922) 

3.0 

2.4 

30.0 

-150.47228(-.48336) 

3.0 

2.5 

30.0 

-150.471 87(-.48308) 

3.0 

2.6 

30.0 

-150.46851(-.47976) 

2.8 

2.4 

30.0 

-150.48517(-.49622)' 

2.8 

2.5 

30.0 

-150.48661(-.49778)' 

2.8 

2.6 

30.0 

-150.48456(-.49578)' 

2.6 

2.4 

30.0 

-150.49880(-. 50978)' 

2.6 

2.5 

30.0 

-150.50155(-. 51265)' 

2.6 

2.6 

30.0 

-150.50031(-. 51147)' 

2.4 

2.4 

30.0 

-150.51004(-.52093)' 

2.4 

2.5 

30.0 

-150.51376(-. 52477)' 

2.4 

2.6 

30.0 

-150.51338(-. 52447)' 

2.2 

2.4 

30.0 

-150.51472(-. 52548)' 

2.2 

2.5 

30.0 

-150.51940(-.53029)' 

2.2 

2.6 

30.0 

-150.51953(-. 53050)' 

2.0 

2.4 

30.0 

-150.50681(-.51743)' 

2.0 

2.5 

30.0 

-150.51256(-.52332)' 

2.0 

2.6 

30.0 

-150.51374(-.52458)' 

1.8 

2.8 

30.0 

-150.48787(-.49866) 

1.8 

2.7 

30.0 

-150.48976(-.50051) 

1.8 

2.6 

30.0 

-150.48948(-.50016)' 

1.8 

2.5 

30.0 

-150.48602(-.49660) 

1.8 

2.4 

30.0 

-150.47799(-.48845) 
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3.2 

30.0 

3.0 

30.0 

2.9 

30.0 

2.8 

30.0 

2.7 

30.0 

2.6 

30.0 

2.5 

30.0 

2.4 

30.0 

2.28 

20.0 

2.28 

20.0 

2.28 

20.0 

2.28 

20.0 

2.3 

20.0 

2.4 

20.0 

2.3 

20.0 

2.4 

20.0 

2.5 

20.0 

2.4 

20.0 

2.5 

20.0 

2.6 

20.0 

2.7 

20.0 

2.8 

20.0 

2.7 

20.0 

2.6 

20.0 

2.8 

20.0 


- 150 . 44722 (-. 45793 ) 

- 150 . 44885 (-. 45959 ) 

- 150 . 44871 (-. 45943 ) 

- 150 . 44747 (-. 45815 ) 

- 150 . 44463 (-. 45524 ) 

- 150 . 43952 (-. 45003) c 

- 150 . 43119 (-. 44160) c 

- 150 . 41833 (-. 42860) c 

- 150 . 45356 (-. 46254 ) 

- 150 . 45171 (-. 46078 ) 

- 150 . 44810 (-. 45737 ) 

- 150 . 44237 (-. 45223 ) 
- 150 . 44238 (-. 45235 ) 
- 150.43891 (-. 44935 ) 

- 150 . 44106 (-. 45183 ) 

- 150 . 44329 (- . 45426 ) 
- 150 . 44167 (-. 45275 ) 

- 150 . 46918 (- . 48023 ) 
- 150 . 47314 (-. 48434 ) 
- 150 . 47284 (-. 48414 ) 
- 150 . 46949 (-. 48084 ) 

- 150 . 48701 (-. 49756) c 
- 150 . 49198 (-. 50255) c 
- 150 . 49487 (-. 50547 ) c 

- 150 . 49399 (-. 50439) c 


2.0 

2.7 

20.0 

-150.49796(-.50834) c 

2.0 

2.6 

20.0 

-150.49940(-.50976) c 

1.8 

2.8 

20.0 

-150.49191(-.50214) c 

1.8 

2.7 

20.0 

-150.49424(-.50444) c 

1.8 

2.6 

20.0 

-150.49390(-.50406) c 

1.6 

3.0 

20.0 

-150.47234(-.48250) d 

1.6 

2.8 

20.0 

-150.47664(-.48671) c 

1.6 

2.6 

20.0 

-150.47239(-.48235) c 

1.4 

3.2 

20.0 

-150.44828(-.45846) d 

1.4 

3.1 

20.0 

-150.44936(-.45948) d 

1.4 

3.0 

20.0 

-150.44932(-.45937) d 

1.2 

3.8 

20.0 

-150.42758(-.43790) d 

1.2 

3.6 

20.0 

-150.42788(-.43830) d 

1.2 

3.4 

20.0 

-150.42671(-.43700) d 

1.2 

3.3 

20.0 

-150.42507(-.43528) d 

1.2 

3.2 

20.0 

-150.42233(-.43247) d 

1.2 

3.1 

20.0 

-150.41812(-.42819) d 

5.0 

2.28 

10.0 

-150.45313(-.46208) 

4.5 

2.28 

10.0 

-150.45063(-.45963) 

4.0 

2.28 

10.0 

-150.44523(-.45435) 

3.5 

2.28 

10.0 

-150.43420(- .44366) 

3.5 

2.3 

10.0 

-150.43401(-.44357) 

3.5 

2.4 

10.0 

-150.42934(-.43941 ) 

3.0 

2.3 

10.0 

-150.41506(-. 42612) 

3.0 

2.4 

10.0 

-150.41441(-.42564) 
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3.0 

2.5 

10.0 

-150.41082(-.42200) 

2.5 

2.3 

10.0 

-150.41874(-.42984) 

2.5 

2.4 

10.0 

-150.43054(-.44196) 

2.5 

2.5 

10.0 

-150.43654(-.44826) 

2.5 

2.6 

10.0 

-150.43810(-.45009) 

2.5 

2.7 

10.0 

-150.43630(-.44849) 

2.2 

2.8 

10.0 

-150.45917(-.46980) c 

2.2 

2.6 

10.0 

-150.46435(-.47497) c 

2.0 

2.8 

10.0 

-150.47161(-.48214) c 

2.0 

2.7 

10.0 

-150.47510(-.48561) c 

2.0 

2.6 

10.0 

-150.47576(-.48625) c 

1.8 

2.8 

10.0 

-150.47780(-.48820) c 

1.8 

2.7 

10.0 

-150.48019(-.49056) c 

1.8 

2.6 

10.0 

-150.47964(-.48997) c 

1.6 

2.8 

10.0 

-150.47433(-.48460) c 

1.6 

2.7 

10.0 

-150.47453(-.48475) c 

1.6 

2.6 

10.0 

-150.47161(-.48178) c 

1.4 

3.0 

10.0 

-150.45726(-.46754) d 

1.4 

2.9 

10.0 

-150.45830(-.46851) <i 

1.4 

2.8 

10.0 

-150.45729(-.46743) <1 

1.4 

2.7 

10.0 

-150.45360(-.46367) <i 

1.2 

3.2 

10.0 

-150.43572(-.44611) <i 

1.2 

3.1 

10.0 

-150.43556(-.44586) d 

1.2 

3.0 

10.0 

-150.43366(-.44387) <i 

1.2 

2.9 

10.0 

-150.42961(-.43972) <i 
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5.0 

2.0 

0.0 

-150.40941(-.41785) 6 

5.0 

2.2 

0.0 

-150.45111(-.45990) 6 

5.0 

2.3 

0.0 

-150.45345(-.46236) fc 

5.0 

2.4 

0.0 

-150.44854(-.45754) 6 

4.0 

2.2 

0.0 

-150.44257(-.45144) 6 

4.0 

2.3 

0.0 

-150.44494(-.45392) 6 

4.0 

2.4 

0.0 

-150.44010(-.44917) k 

3.5 

2.2 

0.0 

-150.43014(-.43916) k 

3.5 

2.3 

0.0 

-150.43272(-.44189) k 

3.5 

2.4 

0.0 

-150.42810(-.43736) 6 

3.0 

2.2 

0.0 

-150.40491(-.41429) fc 

3.0 

2.3 

0.0 

-150.40815(-.41769) 6 

3.0 

2.4 

0.0 

-150.40416(-.41381) 6 

2.5 

2.4 

0.0 

-150.41135(-.42244) fc 

2.5 

2.5 

0.0 

-150.42059(-.43164) fc 

2.5 

2.6 

0.0 

-150.42431(-.43533) 6 

2.5 

2.7 

0.0 

-150.42401(-.43499) 6 

2.15 

2.5 

0.0 

-150.45182(-.46259) fc 

2.15 

2.6 

0.0 

-150.45537(-.46613) fc 

2.15 

2.7 

0.0 

-150.45486(-.46562) 6 

2.0 

2.5 

0.0 

-150.46131(-.47194) 6 

2.0 

2.6 

0.0 

-150.46517(-.47580) fc 

2.0 

2.7 

0.0 

-150.46494(-.47559) i 

1.8 

2.5 

0.0 

-150.46714(-.47757) fc 
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1.8 

2.6 

0.0 

-150.47205(-.48251) 6 

1.8 

2.7 

0.0 

-150.47280(-.48329) 6 

1.8 

2.8 

0.0 

-150.47049(-.48102) 6 

1.6 

2.5 

0.0 

-150.46109(-.47132) 6 

1.6 

2.6 

0.0 

-150.46822(-.47850) 6 

1.6 

2.7 

0.0 

-150.47104(-.48137) 6 

1.6 

2.8 

0.0 

-150.47065(-.48104) 6 

1.4 

3.0 

0.0 

-150.45702(-.46743) d 

1.4 

2.8 

0.0 

-150.45839(-.46863) d 

1.4 

2.6 

0.0 

-150.44870(-.45879) <i 

1.2 

3.2 

0.0 

-150.43782(-.44836) d 

1.2 

3.0 

0.0 

-150.43801(-.44832) i 

1.2 

2.8 

0.0 

-150.42961(-.43972) <i 

1.2 

2.6 

0.0 

-150.40785(-.41778) <i 


° Points run in C 2 * (collinear) symmetry (from II). 

6 Points run in C 20 (edge-on) symmetry (from II.). 
c Points from II. 

d Points run in C, symmetry with (52) active space to be compatable with the 
calculations in II. 
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I. Introduction 


The reaction 


J?2 “1“ O 2 — ► HO 2 "I" S (1) 

is thought to be an important initiation reaction in H 2 combustion. The only 
experimental determination [1] of the rate of reaction (1) derives from an indirect 
measurement of the rate of the reverse reaction. 

HO2 B — + H2 4 " O2 ( 2 ) 

The products of reaction (2) correspond to a triplet surface. However, the reactants 
also correlate with the ground state of H 2 O 2 for the singlet surface. 

HO 2 + H -y H 2 OI -» OH + OH (3) 

Here the product is OH radical. The rate of reaction (2) was inferred from the rate 
of disappearance of HO 2 and the rate of formation of OH. This rate determination 
was carried out at only two temperatures (« 300K and rs 500K). Thus, it is difficult 
to extrapolate the rate to flame temperatures (« 2000K). Current kinetics models 
[2,3] differ significantly in the rate which is used for reaction (2). The uncertainty 
in this rate is sufficient that a theoretical estimate of the rate, even one with some 
severe approximations, would be useful. 

Reaction (2) is exoergic by 58 kcal/mol. In accord with this, the saddle point is 
expected to resemble H + HO 2 and a small barrier is expected to this H abstraction 
process. In the present calculations the force constant matrix at the saddle point 
for reaction (2) is computed from ab-initio calculations and the rate of reaction (2) 
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is computed rising transition state theory (TST) with a one-dimensional tunneling 
correction based on an Eckart barrier [4]. 

The computational method is discussed in Section II, the results are presented in 
Section III, and the conclusions are given in Section IV. 

II. Computational Details. 

The basis sets used in the present study are atomic natural orbital(ANO) basis 
sets [5]. These basis sets are optimal for describing the atomic correlation and 
have very small basis set superposition errors, but are sufficiently flexible to be 
used in molecular calculations at both the SCF and Cl level. The 0 basis set is 
(13s8p6d4f)/[4s3p2dlf] and is described in detail in Ref. 5. The H basis set is 
(8s6p4d)/[3s2pld], and is that developed by Almlof and Taylor [5]. 

Most of the calculations were carried out in C, symmetry. The CASSCF active 
space had 6 active electrons distributed among five a' and one a” orbitals. The 
qualitative character of the CASSCF orbitals is: l-4a' correspond to the 0 Is and 
O 2s levels, which axe not correlated in these calculations. The 5a' orbital is a 2p 
lone pair localized mainly on the far oxygen. 6a' and 7a' are the 00 and OH bond 
orbitals. 8a' is the H Is orbital. 9a' and 10a' are correlating orbitals for the 00 
and OH bonds. The la" orbital is an O 2p lone pair localized mainly on the near 
oxygen and the 2a" orbital is a singly occupied O 2p orbital localized mainly on the 
far oxygen. 

Ten electrons were correlated in the CCI calculations. The reference configura- 
tions for the CCI consisted of following 29 spatial occupations formed from the 5-10 
a' and the 1-2 a" orbitals: 


f 5a n 6a! 7 7a! 7 8o fl 9a 10 10a , ®\ 

5a' 1 6a' a 7a' 2 8a' 2 9a'°10a , ° 
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5a' 1 Qa n Ta n 8a' 1 9a' 1 10a 1 * 
5a' 1 0a' 2 7a! 7 8a' 1 9a'° 10a' 1 



5a n 6a' 1 7a 13 8a' 3 9a 10 10a'7 xfla" a 2a' rl - 


5a' 2 6a' 1 7a*8a' 1 9a' 1 10a'®\ lla" 1 2a" a > 
5a' 3 6a' 1 7a' 3 8a' 1 9a'°10a' 1 
5a' 3 6a' 3 7a' 1 8a' 3 9a' # 10a ,# 

5a* 6a' 3 7a' 1 8a' 1 9a' 1 lOa* 


\5a n Qa n 7a n Sa' 1 9a*lQa! 1 j 


5a*6a*7a'° 8a' 1 9a' 3 10a'°\ 


5a*6a' 1 7a' 1 8a' 1 9a*10a ,# 
5a*6a'°7a*8a' 1 9a*10a ,# I 
5a*6a*7a'°8a' 1 9a' 1 10a*/?< la" 3 2a" 1 
5a*6a' 1 7a' 1 8a' 1 9a' 1 10a rt 


/5a' 3 6a'°7a' 3 8a' 1 9a' 1 10a' 1 1 

^ ^ rt I 


5a*6a*7a'°8a' 1 9a' # 10a* 

5a*6a' 1 7a' 1 8a' 1 9a'°10a* 

5a r3 6a fO 7a' 3 8a' 1 9a'°10a' 3 


/ 


The first group of configurations includes products of single excitations among the 
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a* HO 2 like orbitals and single excitations among the a” orbitals. While the second 
group of configurations includes double excitations among the a' HO 2 like orbitals. 
These groups of configurations were found to be the most important classes of 
configurations for this system for the H + HO 2 and H - HO 2 saddle point regions 
of the surface, and included all the configurations with expansion coefficient greater 
than 0.05 in the CASSCF wavefunction. 

The multireference analog of Davidson’s correction [6] was computed to estimate 
the importance of higher excitations. The correction used in the CCI [7] is AE (1 
— Cq)/Co, where AE is the Cl energy minus the reference energy and Cq is the 
sum of the squares of the coefficients of the reference configurations in the Cl wave 
function. In the CCI program two different estimates of Co are used. The second is 
obtained as defined above as the dot product of the valence portion of the Cl vector 
with itsself , while the first is the dot product of the valence Cl vector with the 
valence portion of the Cl vector. In most cases these two estimates of give very 
similar results and we have normally reported the first estimate. In the present 
case, there is enough difference that we have chosen to use the second estimate, 
which we believe to be more reliable. 

The electronic structure calculations were carried out on the NASA Cray XMP- 
48 and the NAS Cray Y-MP/8-32. All the calculations used the MOLECULE[8]- 
SWEDEN[9] system of programs. The normal mode analysis at the saddle point for 
reaction (2) was carried out using the program SURVIB[10], while the TST theory 
calculations were carried out using code written locally. 

III. Results and Discussion. 

The computed energies are tabulated in the appendix [10]. The computed points 
are adequate to define both the diagonal and off diagonal components of the force 
constant matrix at the saddle point. A polynomial was fit to the computed points 
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using the SURVIB program [10]. This polynomial contained 45 terms through cubic 
in the in plane coordinates ( too , *OHi t- HOO, thh» *uid L HHO) and a quadratic 
and quartic term in the dihedral angle defined by the HOO plane of the H0 2 species 
and the HHO plane formed by the two hydrogens and the near 0. This fit had an 
RMS error < 0.1 kcal/mol. The saddle point geometry, vibrational frequencies, and 
barrier height axe given in Table I. From Table I it is seen that the saddle point 
r oh, too , and L HOO are 0.15 a 0 longer, -0.06 a<> shorter, and 2.7 0 larger than 
the geometry obtained for free HOO [12] using the same basis set and a similar CCI 
calculation. The changes in too and 0 are consistent with the slightly elongated 
OH bond, based on the computed minimum energy path for H atom addition to 
0 2 [12]. The thh on the other hand is « 0.8 a 0 longer than in free OH. The latter 
result is consistent with an early saddle point as expected for a very exothermic H 
abstraction process. The saddle point is coplanar with the far H atom 3.8 from 
being collinear with the OH bond (in a direction away from the far O atom). The 
computed barrier height is 3.63 kcal/mol before correction for zero-point energy. 

Table II shows the potential for varying the angles a and the dihedral angle fi 
between the HOO and HHO planes. From Table II it is seen that the most stable 
structure has the HH bond trans to the OO bond (/3 = 180°), with the 90° dihedral 
angle structure next, and the 0° dihedral angle structure ( ds) highest. Fig. 2 shows 
the energy as a function of the variation of /3 for a fixed at 3.63 °. One interesting 
feature of Fig. 1 is a slight minimum at /3 =0.0°. This feature appears to be real. In 
particular it has been found that it is not an artifact of symmetry breaking effects, 
which are found to be negligable. In a hindered rotor model, the variation in energy 
for a full rotation of the dihedral angle with L HHO fixed at 176.4 0 (a=3.63 ) 
is « 0.07 kcal/mol ( 23 cm -1 ). This energy change may be compared to an out 

of plane bending frequency of « 300 cm -1 if the out of plane motion is treated 
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as a bend. This result reflects the near degeneracy of the bending modes i.e. the 
rotational barrier would be zero if the bend modes were degenerate. None the less 
two models were considered: i) a bend model in which the two bending modes of 
the HH group are nearly degenerate and ii) a hindered rotor model in which one 
bending mode remains, but the other vibration is treated as a hindered rotation. 
These two models lead to somewhat different rates within the transition state theory 
model, as discussed below. 

The vibrational frequencies are also given in Table I. Prom examination of the 
normal mode coordinates, u >2 and o >3 are found to be OH stretch and 00 stretch 
modes, respectively. The frequencies of these modes differ from those reported 
for free HOO [12] by -26 and +6 cm -1 , respectively. This is consistent with the 
early saddle point for reaction (2). uji is found to be an OH stretching mode. The 
frequency for this mode is 1795 cm -1 smaller than for free HO 2 . Examination of the 
normal mode coordinates indicates that this mode involves motion of both H atoms 
of the H-HO 2 complex. Thus, the lower frequency results from a larger reduced mass 
as well as some decrease in the bond force constant due to the stretched OH bond. 
The remaining two real frequencies are 347 and 333 cm -1 and these correspond to 
in plane and out of plane bending motion of the HH group. The near degeneracy 
of these two modes is consistent with the near collinear arrangement of the H-H-0 
atoms at the saddle point H-HO 2 geometry. It should be noted here that the out of 
plane bending mode is significantly anharmonic. The 333 cm -1 frequency results 
from a quadratic-quartic fit of the dihedral angle. A quadratic term only leads 
to a frequency of 887 cm -1 . Finally, the frequency corresponding to the reaction 
coordinate is 1786 i 

Table I also shows the zero-point correction computed within the bend model 
using the HO 2 frequencies reported in Ref. 12. The frequencies in Ref. 12 were 
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obtained with the same basis set and comparable computational level to that used 
in the present work. There may be some inconsistency with the present work in that 
in Ref. 12 only a quadratic polynomial was used as compared to a cubic polynomial 
in the present work. These small effects will not be of great significance here since, 
within the transition state theory model, the reactant vibrational frequencies mainly 
effect the effective barrier, through the zero-point energy, and the computed barrier 
height is probably uncertain by « 1 kcal/mol in the present calculation. Using 
the bend model the zero-point corrected barrier is 2.01 kcal/mol. In the hindered 
rotor model the out of plane vibration is removed, which reduces the zero-point 
energy by half of 333 cm -1 or « 0.50 kcal/mol, leading to an effective barrier of 
~ 1.5 kcal/mol. In the hindered rotor model the moment of inertia for the hindered 
rotation is very small leading to widely spaced rotational levels. Thus, the rigid 
rotor partition function is still used. 

Fig. 2 compares the computed rates for reactions (1) and (2) to the experimental 
data [1] and to rate expressions used in combustion models [2,3]. In the transition 
state theory calculations the H atom and HO 2 each have an electronic degeneracy of 
2 and a spatial degeneracy of 1, while H-HO 2 saddle point structure has an electronic 
degeneracy of 3 and a spatial degeneracy of 1. Computed rate curves are shown for 
the bend and hindered rotation models. In both cases the computed barrier heights 
are used without correction. The effect of tunneling is extimated using an Eckart 
potential. From Fig. 2 it is seen that at low temperatures the rigid rotor model 
gives a larger rate due to the lower effective barrier, but at higher temperatures the 
bend model shows a larger rate. This effect arises from the large contribution at 
high temperatures of the small out of plane bending frequency to the vibrational 
partition function. The bend model is expected to be more reliable, since the in- 
plane and out-of-plane bend frequencies are nearly degenerate - a result which is 
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consistent with the very small a value (i.e. the H-H-0 moiety is nearly collinear). 
Thus, to a first approximation the model is like a linear triatomic. Within this 
viewpoint, the very small barrier to “hindered rotation” shown in Fig. 1 is simply 
indicative of the close proximity of this system to the linear triatomic limit. Given 
this, only the rate predicted with the bend model is considered in the comparisons 
to experiment and current kinetic models. 

Fig. 2 also shows the recommended experimental rate expression. At least at 
low temperatures, where experimental data exists, the experimental rate is larger 
than the theoretical rate. This result casts some doubt on the experimented rate, 
since the assumption of no recrossing and neglect of variational effects in the TST 
calculations should result in overestimating the rate. The rate used at 1000K by 
the Langley model [3] is in good agreement with the TST rate based on the bend 
model. The Lewis kinetic model [2] agrees well with the TST rate based on the bend 
model at high and low temperatures but is below the computed rate at intermediate 
temperatures. Also the TST rate shows a marked curvature which is not present in 
the Lewis model. Taking all this into consideration, the TST rate based on the bend 
model is recommended as the most reliable estimate of the rate for this reaction. 
The rate as a function of temperature computed with this model is given in Table 

III. 

IV. Conclusions. 

Transition state theory calculations including an estimate of tunneling through an 
Eckart barrier and based on an ab-initio potential energy surface have been carried 
out for the rate of the reaction 


B -f BO% — * B 2 Oj 
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This reaction is found to have an early barrier with the H-H-0 moiety very nearly 
collinear. Consistent with this the in-plane and out-of-plane HH bending modes 
are found to be nearly degenerate. TST theory calculations have been carried out 
for this model (bend model) and also for a model in which the bend motion of 
the H-H-0 moiety is treated as a bend and hindered rotation (rotor model). It is 
concluded that the bend model is far more realistic in this system where the two 
components of the bend are nearly degenerate. The rate based on the bend model 
is recommended as the best estimate currently available for the rate of this reaction. 
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Figure Captions. 


Fig. 1. Energy for H-H0 2 as a function of the dihedral angle (/?) with the remaining 
geometrical parameters fixed at the computed H-H0 2 saddle point geometry. 

Fig. 2. Comparison of TST calculations and "experimental” rate data for the 
reaction H 2 + 0 2 — < ► H + H0 2 . 
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Table I. Computed Saddle Point Properties for H 4- H00. 


Geometry 



a 

-3.8 

*HH 

2.238 

roH 

1.991 

TOO 

2.463 

6 

107.1 


Frequencies 


wi 

Ct?2 

U>3 

a>4 

^5 

It? 6 

d 2 ° 

1736 

1391 

1226 

347 

887 

1786 i 

d 2 -d 46 

1736 

1391 

1226 

347 

333 

1786 i 


Barrier Height 
AE k 3.63 

A ZPE -1.62 

AE(corr.) 2.01 

° The dihedral angle is fit with a quadratic term. 
b The dihedral angle is fit with a quadratic and a quartic term. 
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Table II. Computed Energies for H + H00 (a and (3 variation). 


a 


roo 

0 

3.63 

180.0 

2.457 

107.0 

3.63 

135.0 

2.457 

107.0 

3.63 

90.0 

2.457 

107.0 

3.63 

45.0 

2.457 

107.0 

3.63 

0.0 

2.457 

107.0 


12.0 

180.0 

2.457 

107.0 

7.5 

180.0 

2.457 

107.0 

0.0 

0.0 

2.457 

107.0 

7.5 

0.0 

2.457 

107.0 

7.5 

90.0 

2.457 

107.0 


*HH 

roH 

D2 

2.218 

1.997 

-151.039454 

2.218 

1.997 

-151.039415 

2.218 

1.997 

-151.039381 

2.218 

1.997 

-151.039350 

2.218 

1.997 

-151.039362 


2.218 

1.997 

-151.03920 

2.218 

1.997 

-151.03939 

2.218 

1.997 

-151.03945 

2.218 

1.997 

-151.03919 

2.218 

1.997 

-151.03926 
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Table III. Computed Rate for Reaction (2). a 


T 

K 

rate 

300 

4.01 

0.2976E-11 

500 

1.79 

0.6586E-11 

750 

1.32 

0.1355E-10 

1000 

1.17 

0.2340E-10 

1250 

1.10 

0.3625E-10 

1500 

1.07 

0.5206E-10 

2000 

1.04 

0.9216E-10 

2500 

1.02 

0.1427E-09 

3000 

1.01 

0.2026E-09 


° Computed using the bend model (see text). 

j ^ ^ M w w/ ^ ^ | tflv V ( * ^ 1 * / *'' s '* 


15 



Table AI. Computed Energies for H + H00° 


a 

Too 

e 

THH 

TOH 

0.0 

2.52 

104.4 

18.16 

1.84 b 

0.0 

2.52 

104.4 

2.4 

2.2 

0.0 

2.52 

104.4 

2.4 

2.0 

0.0 

2.52 

104.4 

2.4 

1.8 

0.0 

2.52 

104.4 

2.2 

2.2 

0.0 

2.52 

104.4 

2.2 

2.0 

0.0 

2.52 

104.4 

2.2 

1.8 

0.0 

2.52 

104.4 

2.0 

2.2 

0.0 

2.52 

104.4 

2.0 

2.0 

0.0 

2.52 

104.4 

2.0 

1.8 

0.0 

2.52 

104.4 

1.8 

2.2 

0.0 

2.57 

104.4 

2.238 

2.009 

0.0 

2.52 

104.4 

2.238 

2.009 

0.0 

2.47 

104.4 

2.238 

2.009 

0.0 

2.42 

104.4 

2.238 

2.009 

0.0 

2.37 

104.4 

2.238 

2.009 

0.0 

2.468 

100.0 

2.238 

2.009 

0.0 

2.468 

105.0 

2.238 

2.009 

0.0 

2.468 

110.0 

2.238 

2.009 


E(hartree) 

D1 

D2 

-151.03432 

(-.04594) 

(-.04528) 

-151.01734 

(-.03518) 

(-.03194) 

-151.02580 

(-.04070) 

(-.03913) 

-151.02563 

(-.03901) 

(-.03807) 

-151.02038 

(-.04035) 

(-.03579) 

-151.02460 

(-.04060) 

(-.03854) 

-151.02157 

(-.03554) 

(-.03445) 

-151.02464 

(-.04684) 

(-.04081) 

-151.02304 

(-.04035) 

(-.03764) 

-151.01589 

(-.03058) 

(-.02926) 

-151.02874 

(-.05243) 

(-.04544) 

-151.02324 

(-.03897) 

(-.03708) 

-151.02469 

(-.04059) 

(-.03857) 

-151.02526 

(-.04134) 

(-.03917) 

-151.02473 

(-.04103) 

(-.03867) 

-151.02287 

(-.03939) 

(-.03682) 

-151.02388 

(-.03962) 

(-.03767) 

-151.02534 

(-.04147) 

(-.03926) 

-151.02498 

(-.04150) 

(-.03905) 
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0.0 

2.468 

115.0 

2.238 

2.009 

-10.0 

2.468 

106.9 

2.238 

2.009 

-5.0 

2.468 

106.9 

2.238 

2.009 

0.0 

2.468 

106.9 

2.238 

2.009 

5.0 

2.468 

106.9 

2.238 

2.009 

10.0 

2.468 

106.9 

2.238 

2.009 

-2.5 

2.468 

106.9 

2.238 

2.009 

-2.5 

2.37 

117. 

2.238 

2.009 

-2.5 

2.37 

97. 

2.238 

2.009 

-2.5 

2.57 

117. 

2.238 

2.009 

-2.5 

2.57 

97. 

2.238 

2.009 

-2.5 

2.37 

106.9 

2.238 

2.1 

-2.5 

2.37 

106.9 

2.238 

1.9 

-2.5 

2.57 

106.9 

2.238 

2.1 

-2.5 

2.57 

106.9 

2.238 

1.9 

-2.5 

2.468 

117. 

2.238 

2.1 

-2.5 

2.468 

117. 

2.238 

1.9 

-2.5 

2.468 

97. 

2.238 

2.1 

-2.5 

2.468 

97. 

2.238 

1.9 

-2.5 

2.37 

106.9 

2.34 

2.009 


-151.02295 

(-.03989) 

(-.03719) 

-151.02525 

(-.04164) 

(-.03925) 

-151.02539 

(-.04174) 

(-.03938) 

-151.02541 

(-.04169) 

(-.03938) 

-151.02506 

(-.04113) 

(-.03899) 

-151.02541 

(-.04173) 

(-.03940) 

-151.02089 

(-.03840) 

(-.03522) 

-151.01831 

(-.03439) 

(-.03211) 

-151.01835 

(-.03509) 

(-.03258) 

-151.02140 

(-.03666) 

(-.03508) 

-151.02326 

(-.04199) 

(-.03807) 

-151.02277 

(-.03793) 

(-.03610) 

-151.02071 

(-.03807) 

(-.03524) 

-151.02385 

(-.03861) 

(-.03716) 

-151.02067 

(-.03964) 

(-.03570) 

-151.02132 

(-.03684) 

(-.03490) 

-151.02009 

(-.03715) 

(-.03443) 

-151.02264 

(-.03697) 

(-.03576) 

-151.02333 

(-.03931) 

(-.03700) 
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-2.5 

2.37 

106.9 

2.14 

2.009 

-2.5 

2.57 

106.9 

2.34 

2.009 

-2.5 

2.57 

106.9 

2.14 

2.009 


-2.5 

2.468 

117. 

2.34 

2.009 

-2.5 

2.468 

117. 

2.14 

2.009 

-2.5 

2.468 

97. 

2.34 

2.009 

-2.5 

2.468 

97. 

2.14 

2.009 


7.5 

2.37 

106.9 

2.238 

2.009 

7.5 

2.57 

106.9 

2.238 

2.009 

-12.5 

2.37 

106.9 

2.238 

2.009 

-12.5 

2.57 

106.9 

2.238 

2.009 


7.5 

2.468 

117 

2.238 

2.009 

7.5 

2.468 

97 

2.238 

2.009 

-12.5 

2.468 

117 

2.238 

2.009 

-12.5 

2.468 

97 

2.238 

2.009 


7.5 

2.468 

106.9 

2.34 

2.009 

7.5 

2.468 

106.9 

2.14 

2.009 

-12.5 

2.468 

106.9 

2.34 

2.009 

-12.5 

2.468 

106.9 

2.14 

2.009 


7.5 

2.468 

106.9 

2.238 

2.1 

7.5 

2.468 

106.9 

2.238 

1.9 


-151.02352 

(-.04098) 

(-.03789) 

-151.02367 

(-.03906) 

(-.03727) 

-151.02236 

(-.03892) 

(-.03660) 

-151.02183 

(-.03821) 

(-.03578) 

-151.02160 

(-.03952) 

(-.03627) 

-151.02258 

(-.03756) 

(-.03599) 

-151.02163 

(-.03781) 

(-.03567) 

-151.02316 

(-.03970) 

(-.03714) 

-151.02284 

(-.03863) 

(-.03672) 

-151.02307 

(-.03986) 

(-.03711) 

-151.02278 

(-.03884) 

(-.03672) 

-151.02152 

(-.03852) 

(-.03581) 

-151.02174 

(-.03700) 

(-.03538) 

-151.02139 

(-.03856) 

(-.03569) 

-151.02195 

(-.03770) 

(-.03573) 

-151.02556 

(-.04106) 

(-.03918) 

-151.02482 

(-.04163) 

(-.03910) 

-151.02548 

(-.04122) 

(-.03914) 

-151.02477 

(-.04188) 

(-.03912) 

-151.02354 

(-.04125) 

(-.03812) 

-151.02536 

(-.04013) 

(-.03865) 
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-12.5 

2.468 

106.9 

2.238 

2.1 

-151.02352 

(-.04154) 

(-.03818) 

-12.5 

2.468 

106.9 

2.238 

1.9 

-151.02523 

(-.04022) 

(-.03857) 

-2.5 

2.468 

106.9 

1.8 

2.1 

-151.02735 

(-.04882) 

(-.04346) 

-2.5 

2.468 

106.9 

1.8 

2.4 

-151.04332 

(-.06928) 

(-.06103) 

-2.5 

2.468 

106.9 

2.0 

2.009 

-151.02452 

(-.04257) 

(-.03933) 

-2.5 

2.468 

106.9 

2.0 

2.1 

-151.02592 

(-.04602) 

(-.04145) 

-2.5 

2.468 

106.9 

2.0 

2.4 

-151.03369 

(-.06055) 

(-.05145) 

-2.5 

2.468 

106.9 

2.14 

1.9 

-151.02412 

(-.03952) 

(-.03772) 

-2.5 

2.468 

106.9 

2.14 

2.1 

-151.02459 

(-.04342) 

(-.03961) 

-2.5 

2.468 

106.9 

2.14 

2.2 

-151.02411 

(-.04554) 

(-.03998) 

-2.5 

2.468 

106.9 

2.14 

2.4 

-151.02616 

(-.05296) 

(-.04372) 

-2.5 

2.468 

106.9 

2.238 

1.8 

-151.02231 

(-.03628) 

(-.03512) 

-2.5 

2.468 

106.9 

2.238 

2.2 

-151.02192 

(-.04230) 

(-.03742) 

-2.5 

2.468 

106.9 

2.34 

1.8 

-151.02432 

(-.03795) 

(-.03689) 

-2.5 

2.468 

106.9 

2.34 

1.9 

-151.02682 

(-.04126) 

(-.03986) 

-2.5 

2.468 

106.9 

2.34 

2.1 

-151.02315 

(-.04019) 

(-.03739) 

-2.5 

2.468 

106.9 

2.34 

2.2 

-151.01998 

(-.03922) 

(-.03506) 

-12.0 

2.457 

106.99 

2.218 

1.997 

-151.02520 

(-.04161) 

(-.03920) 

-7.5 

2.457 

106.99 

2.218 

1.997 

-151.02539 

(-.04177) 

(-.03939) 

-3.63 

2.457 

106.99 

2.218 

1.997 

-151.02547 

(-.04181) 

(-.03945) 

0.0 

2.457 

106.99 

2.218 

1.997 

-151.02547 

(-.04176) 

(-.03945) 

7.5 

2.457 

106.99 

2.218 

1.997 

-151.02525 

(-.04139) 

(-.03919) 
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° [4s3p2dlf/3s2pld] ANO basis set. CASSCF/CCI calculations correlating 10 elec- 
trons using a (52) CASSCF active space and selected reference configurations (see 
text). 

b This geometry is taken as H 4- H0 2 at infinite separation. 
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